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ABSTRACT 

We present ground-based MDM U-band and Spitzer/IRAC 3.6/im-band photometric 
observations of the 72 representative galaxies of the SAURON Survey. Galaxies in our sample 
probe the elliptical E, lenticular SO and spiral Sa populations in the nearby Universe, both in 
field and cluster environments. We perform aperture photometry to derive homogeneous struc- 
tural quantities. In combination with the SAURON stellar velocity dispersion measured within 
an effective radius (a c ), this allows us to explore the location of our galaxies in the colour- 
magnitude, colour-er , Kormendy, Faber-Jackson and Fundamental Plane scaling relations. 
We investigate the dependence of these relations on our recent kinematical classification of 
early-type galaxies (i.e. Slow/Fast Rotators) and the stellar populations. Slow Rotator and Fast 
Rotator E/SO galaxies do not populate distinct locations in the scaling relations, although Slow 
Rotators display a smaller intrinsic scatter. We find that Sa galaxies deviate from the colour- 
magnitude and colour-cr c relations due to the presence of dust, while the E/SO galaxies define 
tight relations. Surprisingly, extremely young objects do not display the bluest (V — [3.6]) 
colours in our sample, as is usually the case in optical colours. This can be understood in the 
context of the large contribution of TP-AGB stars to the infrared, even for young populations, 
resulting in a very tight (V — [3.6])— <t c relation that in turn allows us to define a strong corre- 
lation between metallicity and a c . Many Sa galaxies appear to follow the Fundamental Plane 
defined by E/SO galaxies. Galaxies that appear offset from the relations correspond mostly to 
objects with extremely young populations, with signs of on-going, extended star formation. 
We correct for this effect in the Fundamental Plane, by replacing luminosity with stellar mass 
using an estimate of the stellar mass-to-light ratio, so that all galaxies are part of a tight, single 
relation. The new estimated coefficients are consistent in both photometric bands and suggest 
that differences in stellar populations account for about half of the observed tilt with respect 
to the virial prediction. After these corrections, the Slow Rotator family shows almost no 
intrinsic scatter around the best-fit Fundamental Plane. The use of a velocity dispersion within 
a small aperture (e.g. i? c /8) in the Fundamental Plane results in an increase of around 15% in 
the intrinsic scatter and an average 10% decrease of the tilt away from the virial relation. 

Key words: galaxies: bulges - galaxies: elliptical and lenticular, cD - galaxies: photometry - 
galaxies: structure - galaxies: stellar content - galaxies: fundamental parameters 



1 INTRODUCTION 

Galaxies are fundamental building blocks of our universe, and our 
knowledge of their distribution, structure and dynamics is closely 
tied to our general understanding of structure growth. So-called 
scaling relations, that is correlations between well-defined and 
easily measurable galaxy properties, have always been central to 
our understanding of nearby galaxies. With high redshift studies 
now routine, scaling relations are more useful than ever, allow- 
ing us to probe the evoluti on of galaxy populations over a large 
range o f lookback t imes (e.g.lBell et alj2004l ; IConselice et alj2005l ; 
Iziegler et al.ll2005l : ISaglia et al.||20icll 

The colour-magnitude relati on (CMR) was already recog- 
nised in the sixties and seventies ide Vaucouleurslll96ll ; ISandagel 
1 1 9721 : 1 Visvanathan & Sandagdl 19771) and has served as an impor- 
tant benchmark for theories of gala xy forma tion and evolution since 
(e.g. lBower etaill 19921: Hell et alJl2004l;lBernardi etai1l2005h . The 
main drivers are thought to be galaxy metallicity, which causes 
more metal rich galaxies to be redder, and age, causing younger 
galaxies to be bluer. Galaxies devoid of star formation are thought 
to populate the r ed sequence, while star-forming galaxies lie in the 
blue cloud (e.g. iBaldrv et al.l 120041) . The dichotomy in the distri- 
bution of galaxies in this relation has opened a very productive 
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avenue of research to unravel the epoch of galaxy assem bly (e.g. 
|Pe Lucia et al]|2004l ; I Andreonll2006l : lArnouts et al J2007l) . 



Since its discovery dDiorgovski & Davisll 19871 : iDressler et al.l 

119871) . the Fundamental Plane (FP) has been one of the most stud- 
ied relations in the literature. Given its tightness, like many other 
scaling relations the FP was quickly envisaged as a distance esti- 
mator as we ll as a correlation to understand how galax i es form and 
evolv e (e.g. ISaglia et aljl 19931: [fargensen et~ai1ll996l: IPahre et al. 



1998; Kelson et al. 2000; Bernardi et al. 2003; van der Wei et al. 



2004 ; JHolden et alj 120051 : iMacArthur et al.ri2009l) . It is widely 
recognised that the FP is a manifestation of the virial theorem for 
self-gravitating systems averaged over space and time with physi- 
cal quantities total mass, velocity dispersion, and gravitational ra- 
dius replaced by the observables mean effective surface brightness 
({jj, e )), effective (half-light) radius (R c ), and stellar velocity dis- 
persion (a). Since velocity dispersion and surface brightness are 
distance-independent quantities, contrary to effective radius, it is 
common to express the FP as log(i? c ) = alog(a) + /3{/i c )+7, to 
separate distance-errors from others. If galaxies were homologous 
with constant total mass-to-light ratios, the FP would be equivalent 
to the virial plane and be infinitely thin, with slopes a = 2 and 
j3 — 0.4 (in the notation used here). By studying the intrinsic scat- 
ter around the FP, one can study how galaxy properties differ within 
the observed sample. 

Some projections of the FP, known e arlier in time, have 
also been widely studied. iKormendvl ( 119771) found that the sur- 
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face brightness density (i.e. mean surface brightness) of a galaxy 
changes as a function of its size. This relation is usually known as 
the Kormendy relation (hereafter KR). The correlation is such that 
larger galaxies have lower surface brightness densities, compared 
to their smaller counterparts. The size-luminosity relation (SLR) is 
widely used to esta blish the size evolution of galaxies as a func - 
tion of redshift (e.g. lTruiillo et alj|2006l ; lvan Dokkum et alj|2008r) . 
Finally, the last projection of the FP that we consider i n this pa- 
per is the Faber- Jackson relation JFaber & Jacksonll 19761 : hereafter 
FJR), which relates the luminosity of a galaxy to its stellar velocity 
dispersion. 

The highest quality and best understood scaling relations in the 
optical/near-infrared are the ones for early-type galaxies. This is be- 
cause observationally they are much simpler than spirals, with less 
complicated star formation histories, and less extinction by d ust, 
and thus tighter scaling relations (e.g. lLaurikainen et alj|2010l) . It 
is for that reason that often the two groups are treated separately 
in physically similar relations: the prime example being the rela- 
tions between the central stellar velocity dispersion and absolute 
luminosity of elliptical galaxies (FJR), and the rotation velo city 
and absolute luminosity of disc galaxies JTullv & Fisheal 19771) . In 
an attempt to unify properties of these two groups, spiral galaxies 
are usually studied in terms of their bulge and disc properties. The 
resemblance of bulges to ellipticals has lead t o their inclusion in 
the scaling relations of early-type systems (e.g . Bender et alJI 19921 : 



different groups. An example of this is provided by the measure- 
ment of nuclear cusp slopes (e.g. Ferrarese et al. 1994; Bvu n et alj 



iKhosroshahi et al.l2000l : lFalc6n-Barroso et al.l2002l) . although they 
often reveal a much larger scatter and show, on average, an offset 
with respect to the relations of early-type galaxies. While this might 
not be surprising due to the mentioned effects of (younger) stellar 
populations and dust, part of the reason might also be the (often far 
from trivial) decoupling of the bulge from the disc. 

These scaling relations exist for galaxy parameters at vari- 
ous radii. While photometric quantities inside one effective radius 
are easy to measure, the inherent limitations of traditional (single- 
aperture or long-slit) spectrographs have restricted the measurement 
of the stellar velocity dispersion to the central regions of galax- 
ies. For instance, papers based on the Sloan Digital Sky Su rvey 
(SDSS) data (e.g. llBernardi et alj|2003l : iGraves & Fabej|201(3h use 
velocity dispersions that have been determined from central galaxy 
apertures, corrected to effective velocity dispersions using stan- 
da rd aperture corrections. In this paper we follow the approach 
of ICappellari et alj (120061 hereafter Paper IV) and make use of 
the panoramic capabilities of SAURON to measure velocity disper- 
sions in circular apertures going out to an effective radius (a B ), and 
present scaling relations for which all the parameters are measured 
within the same aperture. This method offers the interesting possi- 
bility of presenting the spiral Sa galaxies in the same relations as 
early-type E/SO galaxies. When measuring a c from the integrated 
galaxy spectrum galaxy broadening can be caused by intrinsic ve- 
locity dispersions, or by galaxy rotation. Despite this uncertainty, 
<j c will still be a measure of the mass in a galaxy inside R c . In addi- 
tion, these velocity dispersions will not be affected by the presence 
of central discs, which oft en show low velocity dispersions (e.g. 
iFalcon-Barroso et alj|2003r) . Since R e for most of our Sa galaxies 
is much larger than the radius inside which the galaxy bulge domi- 
nates, the scaling relations will give us information about both the 
bulges and the inner discs of the Sa galaxies. This paper tends to 
investigate both issues by combining photometry with integral-field 
spectroscopy for a representative sample of E to Sa galaxies, treated 
in a consistent manner with a homogeneous database and methods. 
The importance of this last point should not be overlooked, as sup- 
posedly standard parameters can vary greatly when measured by 



1 19961 : iGebhardt et atll 19961 : ICarolloet aljl997l ; iRest et alj|200ll) 

With these goals in mind, we have carried out an optical 
spectroscopic survey of 72 representative nearby E/SO galaxies 
and Sa galaxies to one R c , using the custom-designed panoramic 
integral-field spectrograph SA URON mounted o n the William 
Herschel Telescope, La Palma l IBacon et al.ll200lL hereafter Paper 
I). The SAURON representative sample was chosen to populate 
uniform ly Mb- e planes, equa lly divided between cluster and field 
objects jde Zeeuw et alJI 20021 hereafter Paper II). The work in this 
paper builds on previous results of our survey on scaling relations 
in other wavelength domains (Paper IV; |jeong et alj2009l , hereafter 
Paper XIII). The reader is referred to other papers of the SAURON 
survey for results on the ste llar kinematics ( Emselle m et alJl2004D 
and k inematic classification JEmsellem et alj20071 ; ICappellari et alj 



12007 ) of early- t ype g a laxies and their s t ellar populations 
1 Kuntschner et alj 120061 : [Shapiro et all l20ld: iKuntschner et alj 
l2010f) and on the kine mati cs (Falcon-Barro so et alj|2006r) and pop- 
ulation of early spirals ( IPeletier et al .120071) . Hereafter, we will refer 
to them as Paper III, IX, X, VI, XV, XVII, VII and XI respectively. 

We present in this paper homogeneous ground-based V-band 
and Spitzer 3.6/im-band imaging observations of the 24 elliptical 
E, 24 lenticular SO and 24 spiral Sa galaxies of the SAURON rep- 
resentative sample. Aperture photometry (growth curve analysis) 
is carried out to homogeneously derive a number of characteristic 
quantities to which the more complex SAURON integral-field 
observations are compared. We introduce the sample selection, 
biases and completeness in § [2] The observations and basic data 
reduction are presented in §[3] We describe the aperture photometry 
and determination of the spectroscopic quantities in § [4] and § [5] 
Bivariate scaling relations are shown in §|6] while the Fundamental 
Plane relation is specifically addressed in § [7] We summarise 
our results and conclude briefly in § [8] Description of the stellar 
population synthesis models and methods used to derive stellar 
mass-to-light ratios for our galaxies are presented in Appendices lAl 
and|B] Scaling relations showing the dependencies with kinematic 
substructure and environment are shown in Appendix [C] Tables 
with the measured quantities are presented in Appendix IDI 

Throughout the paper we adopt the WMAP (Wilkinson Mi- 
crowave Anisotropy Probe) cosmological parameters for the Hub- 
ble constant, the matter density and the cosmological constant, of 
respectively H — p kms" 1 Mpc" 1 , Q M = 0.24 and Q, L = 0.76 
l lSpergel et al.ll2007T) . although these parameters only have a small 
effect on the physical scales of the galaxies due to their proximity. 



2 SAMPLE SELECTION, BIASES AND COMPLETENESS 

The SAURON sample is designed to be representative of the popu- 
lation of early-type galaxies in the nearby Universe. By construc- 
tion, E, SO and Sa galaxies were selected in equal numbers (24 in 
each group) to populate uniformly the absolute magnitude versus 
ellipticity diagram. Within each morphological class, galaxies were 
chosen to sample the field and cluster environments equally (12 in 
the field and 12 in clusters). The sample therefore consists of 72 
galaxies. The basis for the sam ple selection was th e Lyon/Meudon 
Extragalactic Database (LEDA; |Paturel et all 19971) . 

Besides the astrophysically-motivated criteria, the instrument 
specifications impose further constraints on the sample selection: 
6° < 5 < 64° to avoid instrument flexure; cz < 3000 tans -1 
to ensure that all the lines of interest are in the observed spectral 
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Figure 1. Photometric characterisation of the SAURON survey. Top row: V-band effective radius (R c ), mean effective surface brightness ((fJ. c )) and absolute 
total magnitudes. Bottom row: same quantities in the 3.6/^m band from the Spitz.erRR.AC dataset. The sample is defined su ch that galaxies uniform ly populate 
the desired (B-band) absolute magnitude range and that our integral-field observations reach typically as far as one R e (see lde Zeeuw et al.l2002l) . 



range; Ms < — 18 mag to ensure that all central velocity disper- 
sions (above 75 kms -1 ) can be measured; | 6 |> 15° to avoid 
crowded fields and large Galactic extinctions. These restrictions 
make the SAURON set a representative but incomplete luminosity- 
and volume-limited sample of galaxies (see Paper II). 



In Figure [T] we show some of the main properties of the 
SAURON sample in the V- and 3.6/xm-bands. The photometric 
quantities have been derived as outlined in Section [4] The field- 
of-view (FoV) for a single pointing of the SAURON integral-field 
unit (IFU) is 33 x 41 arcsec 2 , and therefore, as shown in the fig- 
ure, covers up to one effective (half-light) radius R c for most of 
our galaxies. Larger galaxies were usually mosaiced with several 
SAURON pointings to reach the effective radius. As intended in our 
sample selection, our galaxies uniformly populate the desired abso- 
lute magnitude range above the selection cut. 



While Figure \T\ illustrates the limits in size, mean surface 
brightness and luminosity of our sample, it still lacks information 
about potential biases (other than the luminosity) and the level of 
completeness of our sample. Preliminar y checks with larger, more 

ly-typ .. 

after B03; lLa Barbera et alj201(il ; ICappellari et alj201 lh reveal that 



3 OBSERVATIONS AND DATA REDUCTION 



complete samples of early-type galaxies (Bernardi et al 



20031 . here- 



our SAURON representative sample covers rather well the parame- 
ter space defined by these three photometric indicators for (/i c )<21 
mag arcsec^ 2 at V-band. We illustrate this using the g-band cata- 
logue of B03 as a comparison sample in Sj6] and SJ7] This dataset 
consists of around 9000 early- type galaxies up to redshift 0.3 with 
stellar velocity dispersions above 70 km s _1 . 



3.1 MDM dataset 

Part of the imaging data were obtained at the 1.3m McGraw-Hill 
Telescope of the MDM Observatory located on Kitt Peak, Arizona, 
over 5 observing runs totalling 40 nights: 2003 March 25-April 6, 
2003 October 27-November 2, 2004 February 18-25, 2005 April 
11-17 and 2005 November 2-6. The entire SAURON representa- 
tive sample was observed (Paper II). The thin 2048 x 2048 pix 2 
backside-illuminated SITE 'echelle' CCD was used, and an addi- 
tional 128 column virtual overscan region was simultaneously ob- 
tained with every frame. In direct imaging mode (//7.5), this yields 
pixels of 0" 508 x 0"508 and a 17.'4 x 17'A field-of-view, ensuring 
proper sampling of the seeing and plenty of sky around the targets 
for sky subtraction. The seeing was typically 1'.'8 to 2'.' 6 and no ob- 
servation with a seeing above 3'.'0 was used. The readout noise and 
gain were typically 3.0 e~ and 5.7 e~ ADLP 1 . The Hubble Space 
Telescope (HST) F555W and F814W filters were used, similar to 
the Cousins V and / optical filters. Long non-photometric expo- 
sures were obtained during most nights. To reach a sufficient depth 
and allow correction of CCD defects when combining the images, 
our stated goal was to acquire at least three long offset exposures in 
each filter for every object. Exposure times for individual long ex- 
posures were typically 400 s in F555W, although we adjusted both 
the exposure times and the number of exposures based on weather 
conditions. To internally calibrate the photometry, we also acquired 
a short (100 s) exposure of every object during the few truly photo- 
metric nights. 
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3.1.1 Data reduction 

The data reduction of the MDM images was carried out using 
standard procedures in IRAIq The bias was subtracted in two steps. 
First, the average of the overscan region of each frame was sub- 
tracted from each column. Second, a 'master' bias was subtracted 
from every frame. Since the bias was found to be stable, we used 
a run-wide average of (overscan-subtracted) bias frames obtained 
at the beginning and end of each night. Dark current was found 
to be non-negligible and was subtracted using a 'master' dark, 
again resulting from a run-wide average of dark frames obtained 
at night during bad weather conditions. All galaxy and standard 
star exposures were then divided by a 'master' Hatfield frame, 
resulting from an average of twilight frames obtained in each filter. 
Surprisingly, the fiatfields were found to vary from night to night 
by up to 2 per cent, so a night average was used whenever possible. 
This is not a major limitation, however, as we are mostly interested 
in azimuthally-averaged quantities (see §[4]l. All images of a given 
galaxy in each filter were registered using the match routine by 
Michael Richmond (available at http://spiff.rit.edu/match/!, based 
on the algorithm by Valdes et al. ( 1995). Star lists were obtained 
using SExtractor JBertin & Arnoutslll996l) . Position uncertainties 
in the registered images were typically smaller than a tenth of a 
pixel. Individual images were then sky subtracted and combined 
using a sigma clipping algorithm and proper scaling. Independently 
of the photometric calibration, the short exposures are essential 
in the central region of many objects, where long exposures are 
often saturated. The combined images were flux calibrated using 
the photometric transformations determined, night by night, in the 
way explained below. We estimate the limiting surface brightness 
of our V-band data (~25 mag arcsec -2 ) as the surface brightness 
level la above the uncertainties in the determination of the sky. 



3.1.2 F814W images 

During the reduction of the data a few important issues were iden- 
tified in the F814W images. First, the images suffered from fring- 
ing at the 1-2 per cent level. Since no nighttime exposure of blank 
fields was obtained, we attempted to devise an alternative correction 
from standard star exposures, and also from the galaxies' exposures 
themselves. Standard star images, however, proved to have too low 
signal-to-noise ratios S/N, while median-combined galaxy expo- 
sures failed to remove all the galaxy signal, as most of our targets 
cover a substantial fraction of the field-of-view. We therefore could 
not correct for the fringing. Second, and most importantly, stars in 
the F814W long-exposure images showed a faint, but extended halo 
around them. This was mostly noticeable in the saturated stars. This 
problem, commonly termed 'red halo point spread function (PSF)', 
is due to the use of thinned CCDs together wi th other effects re- 
lated to the aging of the telescope coating (see lMichardll2002l for 
an in-depth study on the issue). These cause the PSF in the F 814W 
filter to extend to well over 100 arcsec (see also lWu et alj2005l and 
Ide Jondl2008D . Although one could in principle devise a correction, 
it would be rather uncertain. Since this issue will affect any mea- 
surement with this filter, we deemed the F814W dataset unreliable 
and discarded it from our analysis. 

The effective radii presented in previous papers of the 



1 IRAF is distributed by National Optical Astronomy Observatories, which 
is operated by the Association of Universities for Research in Astronomy, 
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SAURON series (Papers IV, VI, IX, X, iMcDermid et all 120061: 
IScott et alj|2009h were computed from combined F814W HST and 
MDM images before we identified this issue. Nevertheless, the val- 
ues calculated there appear to be 1% smaller than the ones measured 
here in the V-band (thus consistent with the presence of colour gra- 
dients) and with a scatter implying a small error of 8%. Note that 
this issue has no effect on any of the conclusions in those papers and 
has a very small effect on the quantities derived from our integral- 
field data, which have been updated in Paper XV and subsequent 
papers of the SAURON series. 



3.1.3 Flux calibration 

During photometric nights, in addition to galaxies, we also per- 
formed re peated observa tions of stellar fields including standard 
stars from lLandola i ll992l) . covering a range of apparent magnitudes 
and colours; for many fields, the observations were repeated at the 
beginning, in the middle and at the end of the night, in order to 
also calibrate the dependence on airmass. In total, we acquired 84 
stellar fields, with several (3 to ~ 10) standard stars in each one. 
The aim was to calibrate the data using photometric solutions of the 
following form: 

mv,std — WlF555W,ins + ZF555W + fcF555W^4 + (1) 

CF555W (?TlV,std — Wl.std), 

where m st d are the standard magnitudes from lLandoltl J 19921) . mi ns 
is the instrumental magnitude, z the photometric zero-point, k the 
atmospheric extinction coefficient, A the airmass and c the colour 
correction coefficient. In practice, for each of our stellar fields, 
we identi fied the s tandar d stars with the help of the maps pub- 
lished by lLandoltl J1992I) , and for each star measured the mag- 
nitude which enters in equation Q] as mi ns . This was done by 
means of standard IRAF tasks, within the noao.digiphot.daophot 
and noao.digiphot.photcal packages. The sky background was eval- 
uated taking the mode of the intensity in an annular region around 
each star, situated 3-4 times the full- width half maximum (FWHM) 
of the stellar profile away from its peak; the star was sky-subtracted 
and the computed magnitude corrected for aperture effects. All the 
stars affected by scattered light or saturated, and those for which the 
fits performed in IRAF did not converge, were removed. We took 
all the remaining stars with measured instrumental magnitudes and 
solved equation[T]for z, k and c. We then grouped the stars accord- 
ing to the night in which they had been observed and solved again 
equation[T]on a 'per-night' basis, keeping the colour term c -which 
is very close to zero- fixed to the 'global' value determined at the 
previous step and fitting for the airmass term k and the zero point 
z only. Given that the standard magnitudes of the reference Lan- 
dolt stars are in the V^-band (Johnson) filter, the photometry of our 
images is based on that filter. During the flux calibration, we have 
therefore converted our images from the HST F555W to the V^-band 
(Johnson) filter. 

The internal accuracy of our flux calibration is around 0.03 
magnitudes. In order to investigate systematic departures of our 
calibration, we compared our measurements to apparent magni- 
tudes measured with the same aperture on HST/WFPC2 PCI im- 
ag es. The set of galax ies used in the comparison is that published 
by lLauer et al.l 1 120051) and for which sky values are reported. The 
choice of aperture was arbitrarily set to R c /2, with the constraint of 
it being larger than 5" to avoid uncertainties derived from the differ- 
ent PSFs and smaller than 15" to be fully included within the PCI 
chip. The difference between our V^-band magnitudes and those of 
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Figure 2. Contours of the surface brightness of three galaxies in the 
SAURON sample (in 0.25 mag arcsec" 2 steps) in both the V- and 3.6/xm 
bands. The largest isophotes are set to reach the limiting surface brightness 
of the survey in each band (Ri25 and s;21 mag arcsec -2 respectively). The 
orientation of the images is such that north is up and east is left. The Hubble 
morphological type is indicated in the bottom right corner. The iield-of-view 
of the observations is much larger than displayed here. 



Lauer's HST/F555W imaging is better than 0.05 magnitudes rms, 
with a small systematic offset. The MDM flux calibration predicts 
slightly brighter (<0.05 mag) galaxies than HST. This shift is, how- 
ever, within the expected V— F555W zero-point transformation for 
different late-type stellar templates (i.e. it conforms to the dominant 
old population in our galaxies] 2 ]. 



3.2 Spitzer/IRAC 3.6^m dataset 

In an attempt to extend our analysis of the scaling relations to 
longer wavebands than the optical, and to alleviate the loss of the 
F814W MDM observations, we decided to use the homogeneous 
IRAC 3.6^m imaging dataset provided by the Spitzer telescope. 
This dataset has the great advantage of being less sensitive to the 
presence of dust and provides a better tracer for the underlying, 
dominant, predominantly old stellar mass component of galaxies. 

We retrieved the InfraRed Array Camera (IRAC) images of our 
sample galaxies at 3.6/im through the Spitzer Science Center (SSC) 
archive. These archival images cover a significant fraction of the 
SAURON sample, and were acquired in the context of several pro- 
grams. The remaining objects were observed as part of the specific 
program 50630 (PI: G. van der Wolk) during Cycle 5, meant to com- 
plete observations of the SAURON sample in both IRAC and MIPS 
(Multiband Imaging Photometer for Spitzer) bands. We used the 
BCD images, and mosaiced them together using the MOPEX soft- 
ware. This avoided the artificial point sources sometimes present 
in the centre of the PBCD images. Details about the reduction are 
given in the IRAC instrument Handboolo The output mosaics were 
then sky subtracted in the standard way. As indicated in the hand- 
book, a zero-point of 280.9 Jy was assumed for flux calibration of 
the data into the Vega system. Further details regarding the data re- 
duction can be found in Paper XV and in van der Wolk et al. (in 
preparation). 

After the data reduction and flux calibration processes, we es- 
timate the limiting surface brightness of our 3.6/im-band data to be 
~21 mag arcsec -2 . In order to illustrate the quality of our imaging, 
in Figure|2]we show isophotal contours down to the limiting surface 
brightness of a few galaxies in our sample. The consistency in the 
photometric depth of both datasets ensures that our measured pa- 
rameters truly reflect the potential structural changes as a function 
of waveband, and are not affected by poor imaging. 



3.3 Distances 

We have made a comprehensive effort to collect the best available 
distance estimates in the literature for all galaxies in our sample. 
We have assigned distance estimates adopting the following priority 
order in the methods used: 

(i) For 42 galaxies, the distances were obtained with the surface 
brightness fluctuation (SBF ) method by Mei et alj ( 120071) for the 
ACS Virgo Clust er Survey, Cantiello et alj ( 120071) using archival 
ACS imaging, and lTonrv et al.ld200lh using ground-based imaging. 
We subtract 0.06 mag from the latter distance moduli (i.e., a 2.7% 
decrease in distance) to conver t to the same zero poin t as the HST 
Key Project Cepheid distances JFreedman et alfeOOll) . 

(ii) For 1 galaxy (NGC 5308), the distance is derived using 
Supernovae type la luminosities from iReindl et alj 1 120050 ■ sub- 
tracting 0.43 mag from the distance moduli to convert to Ho = 
73kms _1 Mpc -1 . 

(iii) For 10 galaxies within 12° of the Virgo cluster center 
(defined by the galaxy M 8 7) and with helioc entric velocities < 
2500 km s -1 , we followed ICrook et alj d2007f) and assigned the 
Virgo cluster distance modulus of 31.092 mag JMei et alj2007l) and 
an error of 0.082 mag due to the depth of the cluster. 

(iv) For 2 galaxies (NGC 2273 and NGC 5448), the distances 



2 see Table 5.2 in the HSTAVFPC2 handbook for zero-point transforma- 
tions I http://www.stsci.edu/hst/wfpc2 i. 



3 http://ssc.spitzer.caltech.edu/irac/dh/ 
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are ba sed on "Look-Alike" galaxies dPaturell 19841) from lTerrv et alj 
12003), calibrated with the HST Key Project Cepheid distances, 
(v) For 4 galaxies, the distances are based on the correlation be- 



Tully 



1 galaxy lum ii 
et alj j2008h . 



calibrated with the HST Key Project Cepheid 



distances. 

(v i) For 1 galaxy (NG C 5198), the d istance is based on th e D n -a 



tween galaxy lum inosity and linewidth (Tully-Fisher relation) from 
~H), calibn 

laxy (NG C 
(see iDressleret all 19871) relation from lWillick et all Jl997h . adopt- 
ing Ho = 73kms _1 Mpc -1 and an error of 0.40 mag in the dis- 
tance modulus. 

(vii) For the remaining 12 galaxies, the distances are based on 
their observed heliocentric radial velocities giv en by NEEJj, using 
the velocity field model of iMould et alj d200Cf) with the terms for 
the influence of the Virgo Cluster and the Great Attractor. 

The methods (i)-(iii) typically yield errors in the distances of 
< 10%, while methods (iv)-(vi) are expected to be accurate to bet- 
ter than < 20%. Comparing these reliable distance estimates for 
60/72 galaxies in our sample with the distance estimates based on 
the observed redshifts using method (vii), we find a (biweight) dis- 
persion of ~ 23%. Taking into account the average 7% error in the 
accurate distance estimates, we thus adopt for the redshift distances 
a typical error of ~ 22%, i.e. 0.487 mag in the redshift distance 
modulus. Tables lD"T1 and lD2l list the adopted distances as well as the 
sources of the distance estimates. 



4 PHOTOMETRIC QUANTITIES 

One of the main goals of this project is the measurement of ho- 
mogeneous photometric quantities. As in Paper IV, we have opted 
for simple, yet frequently used, methods to derive these values. 
This has the advantage of being well reproducible and of allow- 
ing comparison with a wid e range of studies (e.g. iBurstein et alj 
ll987uJ0rgensen et alj|l992l) . The values measured here are based 
on aperture photometry. For each galaxy, in both the l/-band and 
the 3.6/im-band, we extracted radial pro files on circular ap ertures 
using the mge .fit sectors IDl|j package of lCappel lari (2002). Fore- 
ground stars and nearby objects were masked using the SExtractor 
lists generated for the registration of the images (see § 13. Lit , The 
profiles were flux calibrated and corrected for Galactic extinction 
using the Av and A3. 6,1m values from NED, which are based on 
COBE, I RAS and the Leiden-D wingeloo HI emission maps as dis- 
cussed bv lSchlegel et alj ( 119981) . We have made no attempt to cor- 
rect the observed fluxes or luminosities for internal extinction. 

Throughout the figures of this paper, circles denote E/S0 galax- 
ies and diamonds Sa galaxies. Filled symbols indicate galaxies with 
good distance estimates, whereas open symbols denote those with 
only recession velocity determinations. In blue we highlight Fast 
Rotator galaxies, in red Slow Rotator galaxies (see ^5. It and in 
green Sa galaxies. The special case of NGC 4550, a galaxy known 
to consist of tw o counter-rotating stellar discs of similar mass 
JRix et alJI 1992c Paper X), is marked in yellow. 



4.1 Effective Radii, mean effective surface brightnesses and 
absolute magnitudes 

We have determined the effective radii and mean effective surface 
brightnesses for our sample galaxies by fitting our aperture photom- 



4 http://nedwww.ipac.caltech.edu/ 

5 http://www.ittvis.com/ 



etry profiles with R '" dSersiclll968l) growth curves of the form 

rR 

L(<R) = 2vr / I(R')R'dR' (2) 

Jo 

= 1-nnURl jf^l [2n,b n (R/R c ) 1/n ] , 

with 7 the incomplete gamma function, R c the effective radius, I e 
the effective surface brightness (at R c ), n the Sersi c index describ- 
ing th e curvature of the radial profile, and we adopt dCiotti & Bertinl 
1 19991) 



1 4 46 

b n - 2 n 1 h 

3 405 n 



25515 n 2 



(3) 



which is an approximation to better than 10~ 4 for n > 0.36. 

When fitting the growth curve profiles, generally the inner 
« 10" as well as regions outside lcr above the sky level have been 
ignored. The former avoids potential complications due to the point 
spread function and the latter reduces the uncertainties associated 
with imperfect sky subtraction. We have used the integrated Sersic 
profile only to extrapolate the outermost part of the growth curve to 
infinity and estimate the total galaxy luminosity. After this, we have 
determined R c from the radius where the growth curve profiles are 
equal to half this total luminosity, i.e., equal to L e — L(<R C ). In 
other words, we have not adopted the R a values from the Sersic 
fit, even though they turned out to be very similar to those from 
the growth curve profiles after obtaining L c from the Sersic fit. 
While the surface brightness profiles of galaxies are not perfectly 
described by a R}' 4 law (e.g. lCaon et all 19931 ; iGTaham et aljl9961: 
iMacArthur et alj|2003l) . the growth curves for the vast majority of 
the objects in our sample, typically early-type galaxies, were well 
represented at large radii by n = 4. However, the growth curves 
of galaxies displaying extended discs and intermediate to edge-on 
configurations (i.e. mainly spirals, but also some lenticulars) were 
often poorly described by a de Vaucouleurs law, hence we fitted 
a Sersic law with n < 4 instead. Overall, the adopted Sersic in- 
dices n were the same in both the V-band and the 3.6/im-band. 
The adopted Sersic n values together with all the other photomet- 
ric quantities for the sample galaxies are listed in Tables ID1I and 
ID2I The approach of using an R ' 4 growth curve to extrapolate 
the galaxy l uminosity to infinity was the same adopted by classic 
studie s (e.g. IBurstein et alj 1 19871 ; Ijorgensen. Franx. & Kiaergaara 
1 19951) and by previous papers of our survey on scaling relations 
(e.g. Paper IV). This allows for a direct comparison of our numbers 
with theirs. In addition to the measurements of L c and R c , we have 
computed their uncertainties via MonteCarlo realisations. Besides 
including the uncertainty in the background sky subtraction, we in- 
cluded correlations among the pixels in the images in two different 
ways, providing lower and upper limits to the uncertainties. First, 
we assumed that the dominant source of error is Poisson noise and 
that the pixels are un-correlated when we include the errors in the 
sky, except for scales smaller than the FWHM of the PSF. For those 
scales we defined a correlation length (FWHM/pixel scale) which 
we set to 2.0 after some tests. The choice of this correlation length 
does not significantly change the output uncertainties for values be- 
low 6. Second, we assumed that the pixels are fully correlated and 
that this is significantly higher than the Poisson noise (which is also 
included). We adopted as the uncertainty the one produced by the 
first method throughout the paper. The second estimate is a good 
test to assess the maximum error one could expect in the worse 
possible situation. For reference we show the different error mea- 
surements for R c in Fig[3] 

The total apparent magnitude and corresponding uncertainty 
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Re, [3 .6]/R.v= 0.71 ±0.10 




10 



100 



R eV (arcsec) 



Figure 3. Ratio between i? e ,v an d R a ,[3 .6] f° r our sample of galaxies. 
The solid thick line shows the average ratio (R e [3.6]/-Re,v)- The dashed 
line marks the one-to-one relation. Circles denote E/S0 galaxies, diamonds 
Sa galaxies. Filled symbols indicate galaxies with good distance estimates, 
open symbols those with only recession velocity determinations. In blue we 
highlight Fast Rotators, in red Slow Rotators and in green the Sa galaxies. 
The special case of NGC4550, with two similarly-massive counter-rotating 
disc-like components, is marked in yellow. Black/gray error bars denote the 
minimum/maximum uncertainties in our analysis (see i]4. It 



follows from L c , simply asm = — 2.51og(2L e ). The mean ef- 
fective surface brightness was computed by dividing L c , the to- 
tal luminosity measured within one effective radius, by the area 
of the aperture, A c = 7rR c 2 , and expressing it in magnitudes, 
(l-i.c)~ —2.5 log(L c /^lc). Since the effective luminosity and radius 
are correlated, the uncertainty in the mean effective surface bright- 
ness was derived after first computing the (/i c ) values for each of 

the corresponding Monte Carlo realisations of L c an d R e . 

As already shown by other au t hors ( e. g. IPahrej 1 19991; 



:ad y shown by other au t hors ( e . g. IFahrel I199S1; 
D l2003l ; iMacArthur eTaT] 120031 : lLa Barbera eTakl 



Bernardi et alj 

20041) . R e values in the infrared appear to be, in general, smaller 



than those in optical bands. This is mostly due to the fact that galax- 
ies tend to be bluer in the outer parts and there fore emit less light at 
these longer wavelengths ( IPeletier et al.ll990T) . In Figure[3]we show 
the relation we find between the independently measured R e ,v and 
.Re, [3.6] values; R c ,[3.6] is on average 29% smaller than R c ,v- As a 
results of this, and as shown in §[6]and|7] differences are thus found 
in the scaling relations for the two photometric bands. 



torious difference in our sample is that of NGC 4486 (M87), where 
our measured value in the y-band of 106", contrasts w ith oth er val- 
ues in the li terature (95".lde Vaucouleurs et alj|l99ll; 104".|Pahrel 
Il999l ; 17l". lFerrarese et alj|l994l ; 194". lKormendv et alj2009l) . 

In addition, we show the best-fitting FP relations presented in 
§ 17,11 including data from the same sources. We carry out this com- 
parison to assess whether different methods to estimate the struc- 
tural parameters can affect our relations. In order to minimize un- 
certainties, we have re-derived a c for each source's R c value. In 
addition, if not provided, our estimated distances were used to con- 
vert R e from arcsec to kiloparsecs. In the case of iKormendv et alj 
( 120091) . we do not use their tabulated mean surface brightnesses, but 
instead compute them ourselves from the total luminosity and effec- 
tive radii they provide in their Table 1 (columns 9 and 17). This is 
to mimic as much as possible our procedure to compute the photo- 
metric quantities. The figure shows that even though the structural 
parameters individually might vary among the different sources, in 
combination they agree well within the observed scatter, and thus 
our best-fitting parameters should not be biased in any particular 
way. 



4.3 V - [3.6] colours 

Colour measurements have been widely used in the past to ex- 
tract first order information on the stellar c ontent of galaxies and 
const raint different formation scenarios (e.g. I Whitdl 1 9801 ; ICarlberd 
Il984l) . Here we determine the effective colour 

<y-[3.6]) e = -2.51og(L v (< Rc,v)/L [3 . 6] (< fl., v )) + const. 

(4) 
measured within a R e ,v aperture. The choice of R c ,v was preferred 
over -R e , [3.6] to match the aperture used to extract our SAURON 
spectroscopic quantities (see §[5}. Aperture corrections, as devised 
in the IRAC Instrument Handbook, have been taken into account 
when deriving the colours. We will use the information from these 
parameters, together with absorption line-strength indices, to estab- 
lish the stellar population properties of our galaxies in different re- 
gions of the scaling relations presented here. 

Central colours (e.g. within R c ,v/8) were not derived due to 
the complexity in matching the MDM and Spitzer PSFs. Further- 
more, V — [3.6] colour gradients were also discarded given the 
presence of dust in many of our galaxies, which introduces features 
in the colour profiles that cannot be accounted for with a single lin- 
ear relation. Nevertheless, an in-depth analysis of colour profiles in 
the SAURON sample, using the Spitzer/IRAC 3.6 and 4.5/^m bands, 
will be the subject of study in Peletier et al. (201 1, in preparation). 



4.2 Literature comparison 

In order to test the reliability of our measurements, we have com- 
pared them with published values in the literature. This exercise 
can reveal large differences among sources, mainly depending on 
the depth of the data, photometric band and most importantly 
the methodology used to derive the relevant quantities. In Fig- 
ure [4] we show the compariso n of our estimated values with a 
few references in the literature jd e Vauco uleurs et al.lll99ll ; IPahra 
Il999l ; |jarrett et al.ll200d ; lKormendv et alj|2009l) . When necessary, 
the literature major-axis R c estimates have been transformed to 
geometric-mean radius to compare with our values obtained from 
circular apertures. As shown in the figure, the agreement with the 
different sources is generally good with a typical scatter of about 
0.14 dex in R c and 0.2 mag in apparent magnitude. The most no- 



5 SAURON QUANTITIES 

In addition to the aperture photometry extracted from the MDM 
and Spitzer/IRAC images, we have determined a number of 
quantities from our SAURON integral-field data that are key for the 
analysis presented in this paper. These are parameters describing 
the richness in dynamical substructures and the stellar content of 
the galaxies in our sample. They have been computed following the 
procedures detailed in previous papers of the SAURON survey. Here 
we will briefly summarise the main aspects and refer the reader to 
the relevant papers for a full description of the methods employed. 
For convenience, the final set of spectroscopic quantities is listed 
in Tables[D3]and |D4] 
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All: 0.00 ± 0.20 mag 
P99: -0.02 ± 0.21 mag 
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Figure 4. Comparison of o ur aperture photometry wit h the literature. F our le ftmost panels: _R e and a pparent magnitudes fro m this paper are compared to 
those in the literature (RC3. lde Vaucouleurs et aljl99ll; P99. |Pahrdl 19991 ; J+00, |jarrett et alj2000l ; K+09, lKormendv et al.l2009l) for V- and 3.6/^m-bands. We 
subtracted 0.1 mag to the literature A' -band values (i.e. the mean K— [3.6] colour for old stellar populations) to compensate for the colour difference. Mean 
offset and scatter are indicated for all sources (All) and P99 alone (as it is the only sample available in the two bands for the two quantities). Two rightmost 
panels: the best Fundamental Plane relations derived in § |7.1| showing our data as well as the literature data for our sample of galaxies (see details in Jj |4.2> . 
Dashed lines in the Fundamental Plane relations mark the \u uncertainty. 



5.1 Kinematic quantities 

The stellar kinematics of our sample galaxies have been ex- 
tracted following the pr ocedure outlined in Papers IX and X. 
Briefly, we used pPXF JCappellari & Emselleml 120041) to fit a 
linear combination of stellar templates from the MILES library 
l ISanchez-Blazquez et alj2006i) and derive the best mean velocity V 
and velocity dispersion a for each spectrum in our datacubes. In this 
paper we are mostly interested in extracting the true first two veloc- 
ity moments of the line-of-sight velocity distribution (LOSVD), and 
therefore we deliberately do not fit the higher order moments (/13, 
/14). We use the extracted V and a maps to compute Ar, a param- 
eter that measures the specific angular momentum within i? c ,v and 
that has led to the new kinematical classification of galaxies pre- 
sented in Paper IX. Throughout this paper we will identify as Slow 
Rotators (hereafter SR) those galaxies with Ar< 0.1, and as Fast 
Rotators (hereafter FR) the rest (as all the papers in the SAURON 
series since Paper IX). We note that for this sample, the notation 
is consistent wi t h that based on the improved criterion defined in 
lEmsellem etai] J201lh for the galaxies in the ATLAS 3D sampltH 
JCappellarietalj201lh . 

In order to plot some of the scaling relations in § [6] we have 
measured the mean stellar velocity dispersion within an i? c ,v aper- 
ture. For that purpose, we summed up all the spectra available 
within such a circular aperture and then computed a following the 
same procedure described above. Whenever the aperture was not 
fully contained within the SAURON FoV, we used the velocity dis- 
persion calibration of Paper IV (equation 1) to correct our values. 



6 http://purl.org/atlas3d 
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For some aspects in § 17.51 we have also measured the mean stel- 
lar velocity dispersion within an _R Cj v/8 aperture (cr e ,8). As in other 
papers in the survey (see Papers IV and XVII), we adopt a random 

error of 5% for our measured velocity dispe rsion values. 

Finally, we make use of the results in iKrainovic et al.l ( 120081 
hereafter Paper XII) to describe the level of kinematic substruc- 
ture present in our maps (e.g. inner discs, kinematically decoupled 
cores, kinematic twists). 



5.2 Stellar population quantities 

As well as stellar kinematic quantities, we have also measured line- 
strength indices within -R e ,v- In this paper, in order to be fully 
consistent with the procedures employed to derive the stellar kine- 
matics and to minimize the uncertainties in the absolute calibra- 
tion of the line-strength indices, we have opted to measure the in- 
dices in the recently defined Line Index System (LIS) LIS- 14.0 A 
dVazdekis et"ai]|2010L hereafter VAZ10). This method has the ad- 
vantage of circumventing the use of the so-called Lick fitting func- 
tions for the model predictions, which requires the determination 
of often uncertain offsets to account for differences in the flux cal- 
ibration between models and observations. The only required step 
to bring our flux-calibrated data to the LIS- 14.0 A system is to con- 
volve the aperture spectra to a total FWHM of 14.0 A. The choice 
of LIS- 14.0 A over other proposed systems (e.g. LIS-5.0 A or LIS- 
8.4 A) is imposed by the galaxy with the largest a c in our sample 

(i.e. NGC 4486). 

Besides the standard Lick indices jWorthevlll994l) that can 
be measured within our wavelength range (i.e. H/3, Fe5015, 
Mg6), we have also measured the H/3 index presented in 



10 Falcon-Barroso et al. 
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Figure 5. Relation between the H/3 and H/3 line-strength indices for our 
sample galaxies. Colours and symbols as in FigE] The dashed gray lines 
mark the predictions of the IVazdekis et alj 1201 Of) models for different 
metallicities. ([M/H]=[-2.32, -1.71, -1.31, -0.71, -0.40, 0.0, +0.2]). 
Our sample appeal's to closely follow the predictions at high metallicities. 



ICervantes & VazdekisM2009l) . This new index is similar to the clas- 
sical H/3 index, but it has the advantage of being less sensitive to 
metallicity. For convenience, we show the relation between the two 
indices for our galaxies in Fig. [5] We however warn the reader that 
this relation is necessarily biased by our sample selection. Since 
we are using stellar population models with solar abundance ratios, 
we use the [MgFe50]' indejjjto minimise the effects of a-elements 
over abundances (see Paper XVII). For galaxies with i? c ,v larger 
than the SAURON coverage, we applied the line-strength aperture 
corrections devised in Paper VI. As established in Paper XVII, typ- 
ical random errors for our measured values are 0.1 A, with system- 
atic uncertainties being 0.06, 0.15 and 0.08 A for H/3, Fe5015 and 
Mg6 respectively. We assume the same aperture corrections and er- 
rors as H/3 for the H/3 index. 

Figure[6] shows the location of our galaxies in the H/3 and H/3 
vs [MgFe50]' diagrams in the LIS- 14.0 A system. The figure illus- 
trates the main reasons for adopting the H/3 index over the tradi- 
tional Lick H/3 index: (1) the H/3 is more insensitive to metallicity 
than H/3 down to [M/H]« —0.71, which makes the diagram more 
orthogonal, and (2) the vast majority of our galaxies fall within the 
model predictions, which is crucial for a proper estimation of the 
stellar mass-to-light ratios (7*). 

Throughout this paper we investigate the potential dependen- 
cies on age via the line-strength index H/3 . This is to be able to in- 
clude the Sa galaxies in the same diagrams. While the use of ages, 
metallicities and abundance ratios is in general desired, estimates of 
these parameters from a single-stellar population (SSP) analysis, as 
was done for the 48 E and SO galaxies in Paper XVII, are not rec- 
ommended given the more continuous star-formation activity they 
have experienced (see Paper XI for more details on these and other 
caveats). The use of the H/3 index will instead provide a robust first 
order indication of the presence of young stars in our galaxies. 

As is required in some of the relations we are presenting in 
this paper, we have estimated the 7, of our sample galaxies in both 
the V- and 3.6/im-bands. We dedicate Appendices [A] and [B] to the 
description of the set of models and the methods used to derive 
these values. 




12 3 4 

[MgFe50]' (A) 

Figure 6. Line-strength index relat ions for the SAURON sample galaxies 
in the LIS-14.0A system jVazdekis et alj|20101) . Colours and symbols as 
in Fig. PJ] The line-strength indices presented have been measured within a 
R e v aperture. Grids show model predictions for constant ages (horizontal 
lines with labels in Gyr) and metallicities (vertical lines labelled with [M/H] 
with respect to solar metallicity). The top panel shows the classical H/3 Lick 
index (Worthey 1994) versus [MgFe50]'. The botto m panel shows the same 
diagr am, but with the the recently defined H/3 index ICervantes & Vazdekij 
2009) instead of H/3. 



6 SCALING RELATIONS 

In this section we show some of the classic scaling relations for the 
SAURON sample of early-type galaxies. Although much work has 
been devoted to these relations in the literature, mostly separating 
galaxies by their morphological classification, here we will focus on 
how deviations from scaling relations depend on the kinematic in- 
formation and stellar populations from our integral-field data. Since 
we found no significant correlation with the level of kinematic sub- 
structure or environment in any of the considered scaling relations 
(demonstrated in Appendix|C]l, we focus on differences between the 
SR E/S0, FR E/S0, and Sa galaxies. 

In the following sections we derive linear fits of the form y = 
a x + /3 to all relations, except the Fundamental Plane in § [7] We 
started from the F ITEX^Q routine taken from the IDL Astro-Library 
JLandsmarj 1 1 9930 which fits a straight line to data with errors in 



[MgFe50]'=0.5 X [0.69 X Mgfe+Fe50 1 5] 



8 Based on a similar routine bv lPress et alj 11992b . 
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both directions, which we extended to include possible correlations 
between the errors in both directions. To find the straight line that 
best fits a set of N data points xj and yj, with symmetric errors 
Axj and Ayj and covariance Cov (xj,yj), the routine minimizes 



= E 



(//j 



P? 



A 



obs,j 



+ A? n 



(5) 



where the combined observational error is given by 

A 2 ohS:J = (A Vj f + a 2 {Ax 3 ) 2 -2aCov(x j ,y j ), (6) 

and A; nt is the intrinsic scatter, which is increased until the value 
of x 2 P er degrees of freedom is unity. Next, finding the changes in 
a and f3 needed to increase \ 2 by unity, yields the (1-cr) uncertain- 
ties on the best-fit parameters. The values of Xj are normalized by 
subtracting the corresponding observational quantities per galaxy 
by the median value of all galaxies. This choice for the reference 
value (or pivot point) minimizes the uncertainty in the fitted slope 
a and its correlation with the in tercept j3. The details and benefits 
of this approach are described in lTremaine et al.l ( 120021) ■ 

In deriving the errors Axj and Ayj , we include the uncertain- 
ties in all contributing quantities, i.e., the errors in the distances, 
the photometric quantities, the kinematic quantities (stellar veloc- 
ity dispersion) and the stellar population quantities. Correlations in 
the photometric quantities are taking into account via our Monte 
Carlo realizations of § 14. It ; for example, when deriving the error 
in (ne)= —2.5 log [L<,/(-kR 2 )] , we first compute from all realiza- 
tions of L c and R c the corresponding (fj, c ) values, and then derive 
the error as the (biweight) standard deviation. Similar Monte Carlo 
simulations of all quantities are used to calculate the covariance 
Cov(xj,yj), which can be significant especially when the same 
quantities are used in both Xj and yj , such as R c in the Kormendy 
relation and the distance in the size-Luminosity relation. 

Unless mentioned otherwise, for consistency and hence for 
easy comparison with most publications on scaling relations for 
early-type galaxies, only E/SO galaxies with reliable distance es- 
timates are included in the fits (i.e. 46 objects). The resulting best- 
fit relations are written in the figures themselves, while the best-fit 
parameters and corresponding errors are given in Table Q] The val- 
ues indicated by Afl t represent the scatter around the best-fit rela- 
tion in the vertical direction, i.e., the (biweight) standard deviation 
of yj — a Xj — /3, for the galaxies j used in the fit. The value in 
parentheses is the intrinsic scatter, after subtracting, in quadrature, 
from Afl t the (biweight) mean of the combined observational er- 
rors A bs,j, which, as indicated in equation (§}, include potential 
covariances between the variables. The intrinsic scatter values and 
corresponding error estimates can be found in Table Q] We con- 
firmed that well within these errors, the intrinsic scatter values are 
the same as Ai nt in equation [5] when the value of \ 2 P er degrees of 
freedom is unity. The other quoted scatter values Ar.i are based 
on all galaxies within the family/families indicated by the index in 
square brackets, i.e., also the galaxies with distances based on their 
recession velocities. 



6.1 Relations with colour 

We presented a set of CMRs in different photometric bands in Pa- 
pers XIII and XV. In both cases the most deviant objects were galax- 
ies displaying widespread ongoing star formation. Here we bring a 
new element to these works by investigating the location of galaxies 
with distinct kinematic properties (e.g. SR/FR) in the CMR. 

In Figure [7] (top panels) we plot the CMRs for our sample 



galaxies using the effective colour (V — [3.6]) c . There are two main 
sources that can change the slope and increase the intrinsic scat- 
ter in these CMRs: young stellar populations and dust extinction. 
They have opposite effects. Young populations will, in principle, 
shift galaxies down from the relation (galaxies become bluer), while 
dust will make galaxies appear redder than their stellar populations. 
Our sample contains galaxies which clearly exhibit dust (most no- 
tably Sa galaxies). Since we have made no attempt to correct for 
internal extinction, the position of these galaxies in this CMR can- 
not be used to infer information about their stellar populations. The 
reddest objects (NGC 1056, NGC2273, NGC4220, NGC4235 and 
NGC 5953) in the plotted relations are indeed Sa galaxies with very 
prominent dust lanes. 

The bottom-left panel of Fig. [7j shows the colour versus a c 
relation. Here both the FR and SR families define a much tighter 
sequence than in the CMRs shown in the top panels. Particularly 
striking are the relations defined by the SRs and FRs, with intrin- 
sic scatters of only 0.019 and 0.064 mag respectively, as opposed 
to the Sa galaxies with deviations of about 0.30 mag, similar to 
those of the CMRs. The Sa galaxies systematically populate a re- 
gion above the best-fit relation. Given the relative insensitivity of 
a c to dust (at least much less than absolute magnitudes in the op- 
tical), this diagram is potentially useful to remove dusty objects in 
samples destined for CMR studies. It also confirms that the average 
dust content is larger in spiral galaxies than in earlier-type systems. 
Despite these deviant objects, it is worth noting that there are a few 
Sa galaxies in the region populated by SR/FR galaxies. 

As already shown in Papers XI, XIII and XV, our sample con- 
tains a number of objects with clear signs of widespread star for- 
mation (NGC 3032, NGC 3156, NGC 4150, NGC 3489, NGC 4369, 
NGC 4383, NGC 4405 and NGC 5953). Following the arguments 
outlined above, one might expect these galaxies to display the 
bluest colours (in the absence of dust). Surprisingly, this is not 
the case. In order to understand this behaviour, we have plotted 
in the bottom-right panel of Fig. UJ the effective colour versus the 
H/3 index measured within the same aperture. In a ddition to our 
data p oints, we also include colour predictions from lMarigo et alj 
(2008J, hereafter MAR08]j together with line-strength predictions 
from the recently released MILES model^j (VAZ10). This combi- 
nation of colours and line indices is consistent, in the s ense that both 
predictions are based on the same set of isochr ones JGirardi et alj 
20001) and were computed for a Kroupa ( 120011) initial mass func- 
tion. These isochrones take into account the latest stages of stel- 
lar evolution through the thermally pulsing asymptotic giant branch 
(AGB) regime to the point of complete envelope ejection. These 
models show that the dependence of the (V — [3.6]) c colour with 
age is rather subtle even for young stellar populations, being much 
more sensitive to metallicity than age. This is the reason behind the 
lack of a strong signature of the youngest objects in the CMRs of 
Figure UJ 

Now that we have assessed that metallicity is the main driver of 
the (V — [3.6])o colour, we can deduce from the colour-<r c relation 
that there must be a strong correlation between metallicity and a a . 
We have made an attempt to determine this relation by computing 
the best linear relation between metallicity ([M/H]) and V — [3.6] 
colour in the MAR08 models (V - [3.6] = 0.73[M/H] + 3.28, see 
Appendix [A}, and then substituting that in the colour - <r c relation 
found here. The selected galaxies are predominantly old and thus 



P. 
2( 



9 http://stev.oapd.inaf.it/cmd, version 2.2 

10 http://miles.iac.es 
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Figure 7. Colour-magnitude, colour-<r c , and colour-H/3 relations for the S AURON sample of galaxies. The best-fit linear relations to the E/S0 galaxies with 
good distance determinations are indicated in each panel and are shown with solid black lines. Afl t indicates the observed and intrinsic scatters (the latter in 
parentheses) for the galaxies used to obtain the fit, while Ar 1 represent the measur ed values for t he diff erent families around the best-fit relations. Colours 
and symbols as in F igure [5] The lines in the bottom-right panel correspond to the iMarigo et all J2008T) model colour predictions combined with MILES 
IVazdekis et"al~ll2010l) Lick/IDS line-strength indices for different mefallicities (see §[6j]for details): [M/H]=[— 0.71 (dotted-dashed), -0.40 (dashed), 0.00 
(solid), +0.20 (dotted)]. For reference, some ages (in Gyr) are indicated along the different metallicity tracks. Pea ks seen in the colour tra cks are due to an 
enhancement in the production rate of AGB stars, that causes a transient red phase in the integrated SSP colours (see lGirardi & Bertellill 1998T) - 



are well reproduced by single stellar population models. This yields 
a relation of the form 

[M/H] = 0.88 (±0.07) log(cr e ) - 1.92 (±0.16). (7) 

A similarly strong correlation was presented in Paper XVII 
([Z/H](R e ) = 0. 32 log (er e ) - 0.8 2) from an independent set of 
model predictions dSchiavonl 120070 and methodology (i.e. using 
line-strength indices rather than colours). Note, however, that the 
two relations cannot be directly compared as the one presented in 
Paper XVII is based on non-solar scaled stellar population mod- 
els, while the one derived here is not. In order to compare them 
we have determined the metallicity of our models using the re- 
lation from the MAR08 models and then adjusted it assuming 
[Z/H] = [M/H] ± 0.75[a/Fe], where the factor 0.75 is a con- 
stant that depends on the element partition used for the models 



(R. Schiavon, private communication) and [a/Fe] is taken from 
Paper XVII. The slope of the resulting relation still appears to be 
a factor of 2 steeper than that of Paper XVII. The apparent in- 
consistency between the two determinations might indicate some 
issues in the modelling of the always complex AGB phase. The 
comparison with other relations in the literature (e.g. J0rgensenl 
1 19991: iKuntschneJ l2000l: iTraeer et al l |2000|; iThomas et alJi2005t 
Sanchez-Blazq uez et alj 20061 : 1 Proctor et al.l 120081 : lAllanson et al.l 
20091 : iGraves et alj|2009l) is not straightforward either since, in ad- 



dition to the stellar population models used, they are mostly based 
on central aperture measurements (both metallicity and velocity dis- 
persion). 
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Figure 8. Kormendy (KR) and size-luminosity (SLR) relations for the galaxies in the SAURON sample in the V- and 3.6/xm-bands. Symbols as in Fig. [5] The 
thick solid line is the best fit relation (as indicated in the equation in each panel). Ag t indicates the observed and intrinsic scatters (the latter in parentheses) for 
the galaxies used to obtain t he fit while Ar i represent the measured values for the different families around the best-fit relations. The grey area marks the 
extent of the lBernardi et alj q2003l) sample. 



6.2 Kormendy and size-luminosity relations 

The KR and SLR for the SAURON sample of galaxies are shown 
together in Fig. [8] As already noticed by other authors, the value 
one obtains for these relations critically depends on the adopted se- 
lection criteria. In particular, the Malmquist bias caused by select- 
ing galaxies according to their lu minosities can have a large i mpact 
on the resulting best-fit relation JNigoche-Netro et al.ll2008l) . The 
volume and luminosity range of our survey means we lack small, 
high surface brightness, faint objects (as estimated from a compar- 
ison with the parent sample of galaxies in the ATLAS 3D survey, 
ICappellari et al.l201 lh . Nonetheless, limiting our fit to galaxies with 
distances below 25 Mpc (where our sample does not suffer those 
limitations) results in best-fit relations that are identical within the 
uncertainties. The comparison of our data in the KR and SLR with 
that of B03 shows that the strongest bias introduced in our anal- 
ysis by the sample selection is in mean surface brightness. With 
the exception of NGC 5845 (the galaxy with the smallest R B ), our 
galaxies seem to populate rather homogeneously the area defined 
by the B03 sample for (^o,v)< 21 mag arcsec -2 . This implies a 
shallower KR and steeper SLR relation than in the more complete 
B03 sample. 

The observed scatter around the best-fit relations is consistent 
with that found in others studies for galaxies within the same mag- 
nitude range (e.g. lNigoche-Netro et al.ll2008f) . It appears that there 
is still a significant amount of intrinsic scatter in both relations, with 
the SR family displaying the smallest values. SR galaxies tend to 



populate the high luminosity end (with the exception of NGC 4458), 
while FR extend across the whole luminosity range displayed by 
our sample. In general, Sa galaxies appear to be fainter (by ~0.8 
mag) than SR/FR galaxies of the same size. This result is consistent 
with the observe d strong dependence o f these relations on morpho- 
logical type (e.g. lCourteau et alj2007l) . 

One interesting feature is the lack of large, high surface 
brightness (or small and very luminous) g alaxies. This z one o f 
avoidance (ZOA) was already not ed by Bender et alj 1 19920 
and has been later be tter defined by iD'Onofrio et alj 1 120061) and 
ICappellarietalJ l bOllI) . 



6.3 Faber- Jackson relation 

In Figure|9] we show FJRs for our sample. As in previous relations, 
the SR family tends to occupy the high luminosity end of the 
relations. Sa galaxies, however, deviate from the best-fit relations 
in the sense that they systematically populate high luminosities 
for a given <t (for a c <125 kms" 1 ). Di sc galaxies are often 
characterised by the Tully-Fisher relation JTullv & Fished 1 19771) 
instead of the FJR. This is because in low mass systems the 
maximum rotational velocity is a much better tracer of total mass 
than the traditionally measured central velocity dispersions. The 
use of a large aperture for the velocity dispersion measurement 
(a a ) presented here significantly helps to improve matters by using 
a parameter close to the integrated second moment (see Paper IV). 
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Figure 9. Faber-Jackson relation for the galaxies in the SAURON sample in 
the V- and 3.6^tm-band. Symbols as in Fig. [3] The solid line is the best fit re- 
lation (as indicated in the equation on the top left in each panel). The dashed 
line shows the best fit to the SR family (e.g. predominantly old systems). 
Afj t indicates the observed and intrinsic scatters (the latter in parentheses) 
for the galaxies used to obtain the fit, while Ar .] represent the measured 
values for the different families around the best-fit relations. The grey area 
marks the extent of the lBernardi et all 120030 sample. 



We believe that the observed offset is mainly a stellar population 
effect, whereby Sa galaxies are on average younger (and thus have 
in general lower stellar mass-to-light ratios) than earlier types, 
as already discussed in the literature (e.g. Paper XI). This effect 
is further supported by the best-fitting relation for the SRs (e.g. 
predominantly old systems). As shown in the figure, the slope of 
the relation for the SRs is larger (dashed line) than the original 
fit and in better agreement with B03. Nevertheless, the number of 
SRs in our sample is rather limited and thus complete samples (e.g. 
ATLAS 3D survey) are required to confirm the observed trend. 

The FJR for our sample appears to 'bend' toward low a e val- 
ues at low luminosities. T here is a hint of this f eatur e in the relations 
for the Coma Cluster by iBower et alj i 1993) and |j0rgensen et alj 
(1199a) . but it is somehow not seen in other relations in the lit- 
erature based on muc h larger numbers of galaxies (e.g. B03, 
lLa Barbera et alj|201fj) . This might be because galaxies with low 
velocity dispersions are often removed from the samples. Irrespec- 
tive of selection effects, the slopes derived here are slightly shal- 
lower than previously found (e.g. Pahre 1999). In part this is be- 
cause it is common to use central velocity dispersions, which will 
increasingly deviate from the values reported here (measured in a 
larger aperture) for larger galaxies, as the velocity dispersion gradi- 
ents are steeper in the inner parts of larger galaxies. 



7 FUNDAMENTAL PLANE 

In this paper we adopt the following notation for the FP: 

log(i? e ) = a log(<r„) + /3 (Me) + 7- (8) 

The multivariate nature of the FP makes common least-squares 
minimisation algorithms not suitable to determine the main pa- 
ram eters (a, @, 7). The literature is vast on alternative methods 



Graham & Colles j|l997l ; lLa Barbera e"tafll2000l : ISaglia et alj 



1 200 lh . several of which we tested to yield consistent best-fit pa- 
rameters within the corresponding estimated errors. The FP results 
presented in this paper are determined via an orthogonal fit by min- 
imiz ing the sum of the abso lute residuals perpendicular to the plane 
(e.g. |j0rgensen et al.ll 19960 . As in §|6] we include the uncertainties 
in all contributing quantities and take into account correlations in 
the photometric quantities via our Monte Carlo realizations. 

In addition to the fitting scheme, sample selection biases in 
any of the quantities of Eq. [8] c an have an important impact on 
the resulting coefficients JNigoche-Netro et alj I2009T) . We have 
tested the sensitivity of our best-fit parameters to magnitude and 
distance (the two main selection criteria in our sample). Our 
sample is a good representation, in terms of luminosity, of the 
complete early-type galaxy population up to a distance of 25 Mpc, 
beyond which we lack small, faint, high surface brightness objects. 
Nevertheless, the inclusion of galaxies at larger distances hardly 
alters the best-fit parameters. In terms of luminosity alone, the 
best-fit parameters remain within the uncertainties as long as the 
faintest level in absolute magnitude is brighter than Mv ~ — 19 
and M[3.5] ~ —22 magnitudes. Therefore, while distance has no 
major impact on the best-fit parameters, the rem oval of the fainter 
object s in our sample does. As shown in Fig. 7 of lHvde & Bernardil 
( 120091) . the expected bias in the a parameter due to our sample 
selection in terms of absolute magnitude is below 10%. As in other 
relations we restrict our fit to E/S0 galaxies with good distance 
determinations (i.e. filled SR/FR galaxies in the relations). From 
the comparison of our data with those of B03, we have established 
that the lack of objects with (At e ,v}< 21.0 mag arcsec~ 2 in our 
sample effectively means that we miss galaxies with 7? c > 12 kpc. 
We have checked, using the B03 sample, that this effect does not 
seem to bias our results in any particular way. 



7.1 Classic FP 

In Figure [TO] we present the FP for our sample galaxies in both the 
V-band and 3.6/im-band. Sa galaxies share the same relation de- 
fined by early-type systems and are not clearly displaced below the 
relation as they usually are. We believe that the main reason for 
this behaviour is the way we have derived the quantities involved. 
When including spiral galaxies in the FP, it is customary to measure 
the photometric and sp ectroscopic quantities only within the bulge 
dominated region (e.g. lFalcon-Barroso et alj|2002f) . Here we have 
measured the half-light radius, mean surface brightness and velocity 
dispersion in a consistent and homogeneous way for all galaxies re- 
gardless of morphological type. The importance of also computing 
the velocity dispersion within the same half-light radius was already 
shown in Figure [4] even though the photometric quantities taken 
from different literature sources for galaxies in our sample individ- 
ually might vary substantially, when combined with a consistently 
measured integrated velocity dispersion, they all fall on our best- 
fit classic FP. In §[73] we investigate the effects in the best-fitting 
parameters if non-consistent velocity dispersion measurements are 
used. 
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Figure 10. Edge-on views of the Fundamental Plane relation for the galaxies in the SAURON sample of galaxies in V- and 3.6/jm-bands. Symbols as in Fig. [5] 
The dashed lines mark the la uncertainty (measured in log(iJ c )) considering all objects. The solid line is the best fit relation (as indicated in the equation in 
each panel). Ag t indicates the observed and intrinsic scatters (the latter in parentheses) measured along the log(iJ c ) axis for the galaxie s used to obtain the fit , 
while Ar i represent the measured values for the different families around the best-fit relations. The grey area marks the extent of the lBernardi et alj <2003l) 
sample. 



The best-fit coefficients in both bands are somewhat inconsis- 



tent with some of the most recent works in the li t erature (e.g.lPahre 



1999 



2009 



Bernard! et all 120031; IProctor et al J [20081 : iHvde & Bernardi 
La Barbera et alj 120091) . Although difficult to assess accu 



rately, we believe that the main differences are due to a combination 
of photometric bands employed, fitting methods, and specially 
sample selection (since we include the Sa galaxies, unlike previous 
works). There is nevertheless a large disparity in the literature 
regarding the value of the FP coefficients. There seems to be, 
however, a consensus on the fact that a changes gradually with 
photometric band (i ncreasing with wavelen gth) while /3 remains 
almost constant (e.g. lHvde & Ber nardi 2009). The relations derived 
here share that property. The parameters also deviate from the virial 
theorem predictions (i.e. a — 2 and J3 — 0.4 in the notation used 
here), an effect known as the 'tilt' of the FP. At first sight it might 
seem surprising that the /? coefficient in the F-band, more sensitive 
to young stellar populations than infrared bands, has a similar value 
to that in the 3.6/im-band. It is important to remember, however, 
that the 7 parameter, while assumed constant, is in fact a function 
of the total mass-to-light ratio (AI/L), which is not necessarily 
constant among galaxies. In fact this M/L depends on the stellar 
populations, i.e., one can express the total M/L as a function of the 
stellar mass-to-light ratio 7* and the relative fraction of luminous 
to dark matter. It seems that during the fitting process stellar pop- 
ulation effects conspire to keep j3 constant while affecting a and 7. 



7.2 Outliers and residuals 

Not all the galaxies in our sample appear to follow the main FP rela- 
tion. NGC4382 is a well-known shell galaxy (e.g. lKormendv et alj 
120090 which must have suffered some recent interaction event. 
It contains a rather young nuclei (3.7 Gyr, see Paper XVII), it 
is the younge s t of all the non-CO detected SAURON galaxies 
l lCombes et alj 120071) and it also displays a prominent central 
velocity dispersion dip (see Paper III). NGC 1056, NGC4369, 
NGC4383, NGC 5953 and NGC 7742 are Sa galaxies that contain 
a significant fraction of young stars within R c (see Paper XI). 
NGC 3489 is one of the few early-type galaxies with a strong 



starburst in the last 2 Gyr (see Paper XVII). NGC 7332 was already 
recognise d as an outlier of the FP in the optical and near-IR 
bands by iFalcon-Barroso et alj (2002j). Although it shows signs 
of widespread young populations, it is not at the level of other 
early-type galaxies. Nevertheless, this object is peculiar in th at it 
is rich in kinematic substructure JFalcon-Barroso et al.1120040 and 
might have suffered from a recent interaction with NGC 7339. 
Moreover, the presence of two gaseous counter-rotating discs 
l lPlana & Boulesteixll 19961) suggests that this galaxy might not be 
in dynamical equilibrium. It is interesting to note that the object 
in our sample with the strongest presence of widespread young 
populations, NGC 3032, does not deviate at all from the best fits. 

Individual inspection of the outliers to understand the reasons 
for their unusual location with respect to the main FP strongly sug- 
gests young stellar populations as a common factor. We investigate 
this further by exploring whether the residuals correlate with stellar 
population quantities, such as H/3 shown in Figure QT] A correla- 
tion with the H/3 index indeed exists in both bands. With the ex- 
ception of NGC 3032 noted above, galaxies with young stellar pop- 
ulations are systematically located below the main relation. These 
residuals suggest that, while metallicity gradients may contribute 
to the scatter in the FP, young populations are the dominant factor. 



7.3 Influence of discs 

In order to disentangle the possible effects on the FP due to struc- 
tural non-homology in galaxies from those due to young stellar pop- 
ulations, we have carried out a simple experiment illustrated in Fig- 
ure [12] We have defined a mock early-type galaxy as a spheroid 
with a light profile described by a de Vaucouleurs' law. By con- 
struction, we put this spheroid on a generic FP (black filled circle 
and solid line). In order to assess the impact of an additional disc 
component, we added exponential discs to the spheroid, with scale 
lengths of respectively 0.2, 0.5, 1.0, 1.5 and 2.0 times the R c of the 
spheroid (circles, upward triangles, diamonds, downward triangles 
and squares, respectively). For referen ce, we note that the typical 
values are > 5 for late-type galaxies JBalcells et alj|2007l) . Addi- 
tionally, we set the light of the disc to be 0.2, 0.5, 1.0, 2.0 and 5.0 
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Figure 11. Residuals (along the log(i?, c ) axis) of the Fundamental Planes in 
Fig. I lOl as a function of the H/3 line-strength index. Symbols as in Fig. [3] 



times the light of the spheroid (dark blue, cyan, green, red, and pink, 
respectively). We have analysed the growth curves of all these mock 
galaxies in the same way as the real galaxies presented in this pa- 
per. We have assumed the same <t c for all mock galaxies as its value 
does not seem to depend on the Sersic n used in the growth curves 
of our real galaxies. 

The figure shows that with an increase of the disc light frac- 
tion, the mock galaxies tend to deviate more and more from the 
original FP. This is expected as the increase of light will lead to 
higher surface brightnesses, which in turn will push the objects be- 
low the relation. It is interesting, however, that the deviations only 
seem significant when the light fraction of the disc relative to the 
bulge is above a factor 2, regardless of the size of the disc relative 
to the spheroid. We note the large excursions in the R c and verti- 
cal directions for large disc light fractions. Figure [12] also shows 
the effect of a compact or extended disc relative to the underlying 
spheroid. Compact bright discs tend to shift the galaxies towards 
the left of the relation, while extended ones shift towards the right. 
It is important to realise that the presence of a disc has little influ- 
ence on the end location of a galaxy, unless the disc is rather bright 
relative to the spheroid. 

There are, in principle, other sources for the scatter and tilt 
of the FP that we are not accounting for here in detail: projection 
effects, rotation, structural homology, etc. All those, besides stel- 
lar populations, hav e been studied in detail in the literature (e.g . 
Guzman et alj|l993l: ISaglia et al.l|l993l: [Prugniel & Simienl | l994| : 



Ciotti et al.lll996l: Lforgensen et alJl 19961: [Graham & Colles; 



r3|l994j: 
islll99i 



Prugniel & Simienlll997l: IPahre et al.|[l998l : iMobasher et"al]ll999l : 
Truiillo et alj 120041) , although with different and sometimes con- 



flicting conclusions. Recent studies, however, po int to stellar pop- 
ulations and dark matter as the main drivers JTreu et al.l 120061 : 
ICappellari et al.l2006l : lBolton et al.l20o3 : lGraves & Fabej2010|) . In 
this paper we focus on the effect of stellar populations only. 

Results on the importance of young stellar populations in de- 
termining the location of galaxies in the FP were already high- 
lighted in Paper XIII, based on GALEX observations in the ultra- 
violet regime. Although the analysis in that paper was limited to 
fewer galaxies than the sample presented here, it appears that an 
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Figure 12. Simple experiment to show the effects of adding a disc to a 
spheroid on the location of galaxies in the FP (see § |7.3> . The black dot 
marks the location of a pure spheroid, by construction on the FP. The black 
solid line marks the FP, and for reference dashed and dotted lines denote 
the la and 3<r scatter in the relation shown in Figure [Tol Different sym- 
bols indicate the ratio between the scale length of the added disc (h r ) and 
the effective radius (R c ) of the underlying spheroid. Colours correspond to 
different fractions of disc light (Ldisc) relative to the spheroid (£ S pheroid)- 



important fraction of the tilt and scatter of the UV FPs is due to 
the presence of young stars in preferentially low-mass early-type 
galaxies. Triggered by those results, we have attempted here to take 
a step further and correct the FP for stellar population effects in or- 
der to bring all galaxies into a common FP relation. The results of 
this exercise are shown in the following section. 



7.4 FP corrected for stellar populations 

A first test to assess the importance of young stellar populations on 
the tilt of the FP is to derive the best-fit relation after restricting 
the sample to predominantly old galaxies, i.e. SRs and FRs with 
H/3 < 3.0 A, and good distance determinations (35 objects). This 
simple experiment yields the best-fit parameters a = 1.39, j3 = 
0.35 and 7 = -9.61 in the V^-band and a = 1.63, /3 = 0.34 and 
7 = —8.85 in the 3.6/im-band. The dramatic change compared to 
Fig.[l0]is the sudden increase of the a parameter, while j3 remains 
almost unaltered. This object selection reduces the tilt of the FP 
and brings it much closer to the virial prediction. Nevertheless it is 
worth noting that the two a values (i.e. one for each band) are not 
close to each other, as one might expect if we had fully corrected 
for stellar population effects and hence were left with effects, in 
particular the luminous-to-dark matter fraction, that are insensitive 
to wavelength. 

Considering the residuals in Fig. QT| one can take this exper- 
iment a step further and try to fully compensate for stellar popu- 
lations (young and old) simultaneously. H/3 is not a quantity that 
directly enters the FP equation as defined in Eq. [8] However, as 
stated in Sect. 17.11 the 7 coefficient does depend on stellar popu- 
lations via 7,,. We have made an attempt to estimate the effective 
7* for each galaxy based on the combined MILES+MAR08 mod- 
els, assuming a combination of two single stellar populations as the 
baseline star formation history (see Appendices[A]and[B]l. The rela- 
tions between our computed 7, and observed H/3 index are shown 
in Fig. [T3] one for each photometric band. It appears that H/3 is 
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Figure 13. 7* vs H/3 relations for the SAURON sample galaxies in the V 
and 3.6^tm-band. Best-fitted (log-log) linear relations are indicated in each 
panel and plotted with a dotted line. 7* values have been determined as 
explained in Appendix IS1 



strongly correlated with our measured 7, and thus could be used as 
a rough surrogate for 7* . This behaviour was already shown in Pa- 
per IV for the H/3 index. The fitted relations, despite being strong, 
have two important shortcomings: (1) they over-predict 7, at the 
low H/3 end and (2) they are rather uncertain for H/3 values above 
3.0 A. The first effect is an inherent limitation of the models used, 
i.e. any given set of models predicts a maximum 7, for a given ini- 
tial mass function. The second effect is related to the wide range of 
allowed 7* for a given H/3 value, as parametrised by the set of two 
single stellar population models used here. 

There are a few notable exceptions to the strong correlation 
at 3.6/im. These are objects located very close to the peaks of the 
V — [3.6] colour in the MAR08 models for solar metallicities or 
above (i.e. NGC 1056, NGC4383 and NGC 5953; seeFig.UJ. Their 
H/3 values range from 3.8 to 4.2 A. Since our 7* p.e] predictions 
are made after correcting for the colour difference with j* t v (see 
Appendix^! it is not completely surprising that the 7* J3.6] values 
in that region somehow deviate from the trend defined by neigh- 
bouring objects. However, sudden jumps in 7* of that magnitude 
seem unrealistic for these galaxies, and thus we opted to correct 
their values using the linear 7* vs H/3 relation in that band pre- 
sented in Figure[T3"1 

Considering the nearly linear log-log relations between 7* and 
H/3 , we can re-write the original FP equation in terms of these 
variables: 

log(7? c ) = a'log{a s )+/3'(n e )+'y' + 5'\og(R/3 ), (9) 



log(R c ) = a' log(ff c ) + /3'< M c> +1+5' log( 7 «), 



(10) 



where 7' now depends (mainly) on the luminous-to-dark matter ra- 
tio. 

In Figure[l4]we plot the results of this experiment. In order to 
make a meaningful estimate of the effects of the stellar populations, 
we consider all objects with good distance estimates in the fit (60 
galaxies, as opposed to the fits in Fig. 1 101 where only 46 E and SO 
with good distances were considered). We perform this test in both 



bands. The top row in Fig. [14] shows the resulting fits using the 
classic FP definition (eq.[8]l, whereas the middle and bottom panels 
present the fits including the H/3 and 7* terms (equations [9] and 
UOt . It seems that the inclusion of the Sa galaxies in the fits of the 
original FP relation (top panels) results in smaller slopes, which 
translates into larger departures from the virial expectations. 

Our first attempt to bring deviant objects closer to the relation 
is to add the H/3 term in the FP relation. Unfortunately, leaving 
all the input parameters unconstrained results in very poor fits that, 
while successful in bringing back the most deviant galaxies, display 
much larger scatters than previous relations derived in this paper 
and, more importantly, have unreasonable coefficients. In order to 
get meaningful fits, we have instead fixed the H/3 terms to the co- 
efficients expected from the fits to the original FP residuals versus 
H/3 shown in Fig. QT] The results of fixing this term are shown 
in the middle panels. These fits are able to bring the young objects 
much closer to the relations and at the same time to slightly reduce 
the intrinsic scatters. The a 1 coefficients have increased with re- 
spect to the original relation, getting them a step closer to the virial 
prediction, while the (3' values have only changed within the un- 
certainties (~ 0.01). The main disadvantage of H/3 as a substitute 
to 7» can be seen in the case of NGC 3032. As already shown in 
Fig. [T3] this galaxy is a deviant object in the 7* vs H/? relation, 
and this translates in it being pushed out of the relation (its mea- 
sured H/? value over-predicts its true 7,). 

Finally, in the bottom panels of Figure [14] we add the 7, term 
in the FP relation. In this case the unconstrained fits result in phys- 
ically meaningful coefficients, while at the same time bringing the 
most deviant points closer to the relation. The best-fit parameters 
reduce the FP tilt compared to the original relations by almost 50% 
in each band. It is important to remark that best-fit coefficients in 
both bands are now consistent within the errors, as expected if we 
had corrected for the effects of stellar populations. The coefficients, 
at least in the 3.6^im-band, are also within the uncertainties of those 
obtained by fitting the old galaxies only (see above). While in the 
V-band most of our young objects get closer to the relation, in the 
3.6/im-band this is not the case (even though we have corrected 
their 7* estimates based on Fig. [13] as explained above). In fact, 
these objects cannot be brought back to the 3.6/im-band relation 
even with the most extreme 7* values allowed by our stellar popu- 
lation fitting procedure. 

As an additional exercise (not shown here) we have studied 
the effect of using the (V — [3.6]) e colour, instead of H/3 or 7,, in 
reducing the tilt of the FP. As expected, given the poor sensitivity 
of this colour to age (as shown in Fig. [7j, this quantity does not 
help reducing the tilt. At the same time it also demonstrates that 
metallicity cannot be a major player in producing the tilt (since this 
colour correlates strongly with metallicity, see Appendix [At. 

It is important to remind the reader that we have only cor- 
rected for the effects produced by stellar populations, and partly 
by rotation (by using a c ). Still there are a number of other factors 
that influence the final location of galaxies in the FP. The complex 
manner in which star formation and dust compensate each other 
might be at the heart of the location discrepancies of some of the 
objects in the V- and 3.6/xm-band FPs. Nevertheless, the results 
of our exercise demonstrate that it is possible to correct for stellar 
population effects and bring most galaxies to a common relation, 
regardless of their morphological type and photometric band em- 
ployed by including sensible estimates of the s tellar mass-to-light 
ratios. This approach was previously adopted bv lPrugniel & Simienl 
( 119961) in a sample of elliptical galaxies and has been more recently 
exploited by other groups on much larger samples, though still 
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Figure 14. Edge-on views of the Fundamental Plane relations for the galaxies in the SAURON sample in the V- and 3.6/^m-band. Symbols as in Fig. [5] (Top 
row) Original FP relation. (Middle row) Original FP relation with a H/3 term. (Bottom row) Original FP relation with a 7* term. The dashed lines mark the 
la uncertainty (measured along the log(i? c ) axis) considering all objects. The solid line is the best fit relation (as indicated in the equation in each panel). 
Afl t indicates the observed and intrinsic scatters (the latter in parentheses) measured along the \og(R c ) axis for the galaxies used to obtain the fit, while Arj 
represent the measured values for the different families around the best-fit relations. 



restri c ting the analysis to ear l y-type systems (e.g. Allanson et al.l 
120091 ; iHvde & Bernardill2009l ; braves & FaberJl2010h . Our results 
are largely consistent with those. 

7.5 Scatter in the FP 

One of the most striking features observed in our fits of the FP is 
the very tight relation defined by the SR galaxies in both bands. 
This is not a totally unexpected result as SR galaxies are uniformly 
old, but it is remarkable how the trend is kept even for smaller and 
fainter galaxies (as SR galaxies extend over the whole range in 
luminosity of our sample). The FR family displays slightly larger 
rms (A[ PR j) values. The scatter of the Sa galaxies appears to be 
the largest of the cases we have studied here. The different panels 
in Fig. [TO] and [14] show the observed as well as the intrinsic scatter 
(within brackets) for each family. For easily comparison with other 
works in the literature we choose to provide them measured along 
the log(i? ) direction. 

The first thing to notice for the classic fits in Fig.|10|andll4lis 



the differences between the values reported. The scatters for the Sa 
galaxies in the first figure are larger than in the second. Conversely 
the rms for the SR/FR families is smaller in the first figure than in 
the second one. This is simply due to the fact that for determining 
the FP we fit only the SR/FR galaxies in Fig. [TO] while Fig.ll4lin- 
cludes the Sa galaxies as well. If we focus on this last figure we 
observe, apart from a change in the slopes, an improvement (i.e. a 
decrease) of the scatter for all families when we include the H/3 
and 7* terms to correct for the effects of young populations. It is 
also interesting to note that once corrected for this effect, the intrin- 
sic values found for the different families indicate that while FR and 
Sa galaxies are consistent within the uncertainties, they are clearly 
different from those of the SRs. This finding emphasises the re- 
sults presented in Paper IX and subsequent papers in the SAURON 
series suggesting, based on our kinematic classification, that these 
may be intrinsically a different kind of galaxies. This clear distinc- 
tion is not found between morphologically classified elliptical and 
lenticular galaxies. 

The comparison of the scatter presented here with others in 
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Figure 15. Edge-on views of the Fundamental Plane relations for the galaxies in the SAURON sample in the V- and 3.6/^m-band. Here <r c /8 (cr ,8) is used 
instead of <r e . Symbols as in Fig. [5] (Top row) Original FP relation. (Middle row) Original FP relation with a H/3 term. (Bottom row) Original FP relation 
with a 7„ term. The dashed lines mark the la uncertainty (measured in log(R, c )) considering all objects. The solid line is the best fit relation (as indicated in 
the equation in each panel). Ag t indicates the observed and intrinsic scatters (the latter in parentheses) measured along the log(ii c ) axis for the galaxies used 
to obtain the fit, while Ar .] represent the measured values for the different families around the best-fit relations. 



the literature has to be treated with caution as, while most works 
concentrate in E/S0 galaxies, the different samples were selected 
with very different criteria in terms of luminosity or velocity dis- 
persion. In addition, the methods used to fit the relation can af- 
fect this comparison, although the minimisation of the residuals 
via an orthogonal fit is the most widely adopted. There is a gen- 
eral agreement that the observed scatter is consistent among pho- 
tometric bands, a result we share. The typical observed scatter 
for a general population of E/SO galaxies is around 0.10 dex in 
log(iL) (e.g.|j0rgensen et al II 19961 ; IPahre et al.|[l99& IZibetti et all 
|2002| ; lLa Barbera et al.ll20100 . with decreasing values (~0.07 dex) 
as soon as low velocity dispersion (a < 70km s _1 ) or fainter 
galaxies (M r > —23 mag) are removed from their samples (e.g . 



J0rgensen et al .ill 9961 : iBemardi et al.l r2003; La B arbera etaHl2008l : 



Gargiulo et alj|20091) . The values measured here for the SR/FR are 
in good agreement with those in the lit erature for gene r al sam ples 
of E/S0 galaxies. In the specific case of |j0rgensen et al.l |l996J), the 
lowest observed scatter, 0.047 dex, is measured for galaxies with 



no presence of discs and M r < —23.1 mag. They report that this 
value is consistent with the measurement errors. The closest set of 
galaxies in our sample matching those constraints is the SR fam- 
ily. For this group, we find a very small amount of intrinsic scatter 
left, 0.021 ± 0.029 dex and 0.012 ± 0.050 dex in the V-band and 
3.6mu-band respectively. The uncertainties in the our estimates sug- 
gest that the scatter around the best-fit FP are fully explained by the 
observational errors. 

Another interesting point to discuss is whether the use of a 
large aperture velocity dispersion, instead of a central value, helps 
to reduce the intrinsic scatter in th e rel ation. This was nicely dis- 
cussed in lGraham & CollessJ ( fl997h and lBusarello et alj l fl997l) . us- 
ing samples of elliptical and lenticular galaxies and dynamical mod- 
els with some assumptions on the shape of the velocity ellipsoid. 
The SAURON spectrograph gives us now the chance to repeat this 
exercise using real data. We have measured the velocity dispersion 
within i? c .v/8 (cr ei s) for our sample of galaxies and then fit the FP 
in the same manner done in Fig. [14] This is presented in Fig. Q3] 
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The different panels show that the use of a small aperture veloc- 
ity dispersion affects the resulting fits in two important ways: (1) 
increasing the intrinsic scatter of the relations for all families of 
objects, (2) the measured parameter a has lower values than those 
obtained when a a is used. The first point is mainly due to the fact 
that with the smaller aperture we are probing regions of the galaxies 
where intense star formation is taking place or inner discs exist, and 
thus the measured velocity dispersion is more representative of ar- 
eas around the nuclei than the galaxy as a whole. The second effect 
is likely produced by the fact that a c .s systematically departs from 
CT e (i.e. it is larger) for increasing er values. This, together with 
the expected relation between R c and a Q (given the SLR and FJR), 
helps to reduce the measured slope of the FP relation. Overall, we 
estimate an increase of around 15% in the intrinsic scatter and an 
average 10% decrease in a as a result of using a central value for 
the velocity dispersion term. It is also worth noting that the scat- 
ter of the Sa galaxies has not increased as much as for the SR and 
FR families. This supports the idea that star formation in early-type 
systems takes place predominantly in their central regions, while 
for Sa galaxies this can happen over a much larger region (see also 
Paper XI), and therefore it does not necessarily make much differ- 
ence to measure the velocity dispersion in a small or large aperture. 
Alternatively, this result could be due to the shallower nature of the 
stellar velocity dispersion profiles of the Sa galaxies compared to 
those of the SR/FR galaxies. 



8 CONCLUSIONS 

In this paper we report results from photometric follow-up con- 
ducted in the context of the SAURON project. We use ground-based 
MDM V^-band and Spitzer/IRAC 3.6/im-band imaging to charac- 
terise our sample of E, SO and Sa galaxies. We perform aperture 
photometry to derive homogeneous half-light radii, mean effec- 
tive surface brightnesses and total magnitudes. Combined with the 
SAURON integral-field spectroscopic observations, this allows us to 
explore and understand the location of the galaxies in the main scal- 
ing relations as a function of the level of rotation, kinematic sub- 
structure, stellar populations and environment. 

A number of conclusions can be derived from the work pre- 
sented and are summarised in the following points: 

• The level of kinematic substructure (as determined from our 
stellar kinematic maps) or the environment do not lead to a pre- 
ferred location in any of the scaling relations investigated. How- 
ever, this is subject to the potential biases introduced by our sample 
selection. Our cluster environment is defined by the Virgo Cluster 
and Leo group. 

• The Slow Rotator (SR) galaxies define tighter relations than 
Fast Rotator (FR) galaxies. While this is not totally unexpected (i.e. 
SR are uniformly old systems), it is remarkable how the trend is 
kept both for small and large R a . We find that the SR and FR galax- 
ies do not populate distinct locations in the scaling relations. The 
Sa family is the main contributor to the scatters. 

• Sa galaxies deviate from the colour-magnitude and colour-<7 c 
relations due to the presence of dust. SR/FR galaxies instead define 
very tight relations. 

• Surprisingly, extremely young objects do not display the bluest 
(V — [3.6]) colours of our sample, as is usually the case for op- 
tical colours. This can be understood in the context of the large 
contribution of TP-AGB stars to the infrared even for young pop- 
ulations, which makes the (V — [3.6]) colour almost insensitive to 
age for populations above ~ 1 Gyr. This effect results in a very tight 



( V — [3.6] ) — ct c relation that allows us to define a strong correlation 
between metallicity and a B . 

• A large number of Sa galaxies appear to follow the main rela- 
tions defined by earlier-type systems, and are not necessarily offset 
as previously found. We believe this is due to the use of a e (the stel- 
lar velocity dispersion integrated within R c ), which is not as much 
influenced by the presence of young inner discs and accounts to a 
large extent for rotation in the galaxies. 

• FR and Sa galaxies that appear offset from the relations cor- 
respond mostly to objects with extremely young populations and 
signs of on-going, extended star formation (as already highlighted 
in Papers XIII and XV). For the specific case of the Fundamen- 
tal Plane, we made an attempt to correct for this effect so that all 
galaxies are part of a tight, single relation. Once this is done, the FP 
coefficients in the V- and 3.6/im-bands are the same within the un- 
certainties. The new estimated coefficients suggest that differences 
in stellar populations account for about 50% of the observed tilt 
with respect to the virial prediction. 

• The observed scatter of the SR family around the Fundamental 
Plane is smaller than that of the FR or Sa galaxies. After correct- 
ing for stellar populations, the SR family shows almost no intrinsic 
scatter around the best-fit Fundamental Plane. 

• The use of a velocity dispersion within a small aperture (e.g. 
-R c /8) in the Fundamental Plane results in an increase of around 
15% in the intrinsic scatter and an average 10% decrease of the tilt 
away from the virial relation (as a result of the correlation between 
R e and cr e , given by the size-luminosity and Faber- Jackson 
relations). 

We have deliberately focused on the role of stellar popula- 
tions on the main scaling relations of early-type galaxies. The ef- 
fects other galaxy properties exert on them (e.g. non-homology, 
anisotropy and dark matter variations) will be studied in ongoing 
and upcoming integral-field spectroscopic surveys. 
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Table 1. Scaling relations best-fit parameters 



to 



Relation 


Band 


a 





7 


8 


A fit 


A[SR] 


A [PR] 


A[ S a] 


N 


Sample 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


Colour-magnitude 


V 


-0.06 ±0.02 


3.30 ±0.02 


- 


- 


0.092 ±0.010 


0.067 ±0.009 


0.101 ±0.011 


0.245 ±0.010 


46 


E/S0 


On 

S 


Colour-cr c 


- 


0.64 ±0.07 


3.27 ±0.01 


- 


- 


0.057 ±0.011 


0.019 ±0.017 


0.064 ±0.010 


0.295 ±0.010 


46 


E/S0 


Co 


Kormendy 


V 


2.05 ±0.22 


19.75 ±0.07 


- 


- 


0.447 ± 0.061 


0.400 ± 0.029 


0.464 ± 0.072 


0.743 ±0.079 


46 


E/S0 




Size-luminosity 


V 


-2.84 ±0.22 


-20.73 ± 0.07 


- 


- 


0.464 ± 0.028 


0.461 ±0.026 


0.489 ±0.028 


0.807 ±0.040 


46 


E/S0 


a 


Faber-Jackson 


V 


-4.99 ±0.71 


-20.61 ±0.11 


- 


- 


0.662 ± 0.071 


0.672 ±0.057 


0.644 ± 0.088 


0.811 ±0.193 


46 


E/S0 





Fund. Plane 


V 


1.15 ±0.05 


0.34 ±0.01 


-8.81 ±0.19 


- 


0.067 ±0.015 


0.044 ± 0.016 


0.075 ± 0.017 


0.161 ±0.037 


46 


E/S0 


!^ 


Fund. Plane 


V 


1.08 ±0.04 


0.33 ±0.01 


-8.56 ±0.18 


- 


0.078 ±0.018 


0.051 ±0.015 


0.070 ±0.018 


0.149 ±0.038 


60 


E/SO/Sa 


^ 


Fund. Plane + Hftj 


V 


1.30 ±0.04 


0.32 ± 0.01 


-9.22 ±0.19 


0.77 


0.073 ±0.016 


0.043 ± 0.014 


0.077 ±0.016 


0.116 ±0.038 


60 


E/SO/Sa 


; — 1 


Fund. Plane + 7* 


V 


1.56 ±0.13 


0.32 ±0.01 


-9.03 ±0.31 


-0.55 ±0.13 


0.064 ± 0.031 


0.021 ±0.029 


0.070 ± 0.029 


0.088 ±0.051 


60 


E/SO/Sa 




Colour-magnitude 


[3.6] 


-0.06 ±0.01 


3.29 ±0.01 


- 


- 


0.087 ±0.010 


0.060 ±0.010 


0.095 ±0.011 


0.234 ±0.010 


46 


E/S0 




Kormendy 


[3.6] 


1.98 ±0.27 


15.94 ±0.08 


- 


- 


0.519 ±0.082 


0.464 ± 0.042 


0.541 ± 0.085 


0.745 ±0.097 


46 


E/S0 




Size-luminosity 


[3.6] 


-2.93 ±0.26 


-23.79 ± 0.08 


- 


- 


0.541 ±0.029 


0.396 ± 0.024 


0.553 ±0.120 


0.659 ±0.128 


46 


E/S0 




Faber-Jackson 


[3.6] 


-5.62 ±0.69 


-23.69 ± 0.10 


- 


- 


0.640 ± 0.069 


0.669 ± 0.054 


0.618 ±0.085 


0.911 ±0.187 


46 


E/S0 




Fund. Plane 


[3.6] 


1.37 ±0.05 


0.33 ±0.01 


-8.06 ±0.19 


- 


0.071 ±0.014 


0.037 ±0.012 


0.083 ±0.017 


0.224 ±0.035 


46 


E/S0 




Fund. Plane 


[3.6] 


1.23 ±0.04 


0.32 ±0.01 


-7.49 ±0.16 


- 


0.099 ±0.016 


0.053 ±0.012 


0.084 ±0.017 


0.191 ±0.036 


60 


E/SO/Sa 




Fund. Plane + H/3 


[3.6] 


1.47 ±0.04 


0.31 ±0.01 


-8.35 ±0.18 


0.96 


0.082 ± 0.014 


0.044 ±0.011 


0.080 ±0.015 


0.132 ±0.036 


60 


E/SO/Sa 




Fund. Plane + 7* 


[3.6] 


1.61 ±0.12 


0.32 ± 0.01 


-8.39 ±0.37 


-0.60 ±0.15 


0.057 ±0.041 


0.012 ±0.050 


0.065 ± 0.036 


0.096 ± 0.057 


60 


E/SO/Sa 
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NOTES: 

(1) We use the following notation for each of the scaling relations: 

Colour-magnitude: (V - [3.6]) e = a(M+24) + for the V-band and (V - [3.6]) e = a(M+21) + /3 for the 3.6/xm-band 

Colour-cr c : (V — [3.6]) c = alog(o- c /<7 rc f) + /3, with cr rc { = 150. 

Kormendy relation: {/i e )= ol log(i? c /iJ rc f ) + /3, with iJ rc f equal to 2.4 and 1.7 for the V- and 3.6/im-band respectively. 

Size-luminosity: M = a log(R c / R rc { ) + /3, with R IC { equal to 2.4 and 1.7 for the V- and 3.6/im-band respectively. 

Faber-Jackson: M = a log(cr /<T rc f ) + 0, with tr rc f = 150. 

Fundamental Plane: log(i? e ) = ol log(cr ) + /3{/i e }±7 

Fundamental Plane + H/3 Q : log(-Ro) = a log(cr c ) + /3(^ c )+7 + <51og(H/3 ) 

Fundamental Plane + 7*: log(J? c ) = alog(cr ) + /3(/i c )+7 + (51og(7«) 

Note that the <5 parameter for the Fundamental Plane + H/3 term has been held fixed (see S\7A\ . 

(2) Photometric band of the fit. (3-6) Best-fit parameters and uncertainties. 

(7-10) Intrinsic scatter and uncertainties for the fit and SR/FR/Sa families. Values for the colour-magnitude, colour-cr , size-luminosity and Faber-Jackson are expressed in magnitudes, while in mag 

arcsec -2 for the Kormendy relation. For the Fundamental Plane, the intrinsic scatter is given in dex along the log(i? c ) direction. 

(1 1-12) Number of objects in the fit and type of galaxies fitted. Only galaxies with good distance estimations have been included in the fit. 
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APPENDIX A: STELLAR POPULATION MODELS 



As described in Sect, |6.1| we use a combination of the new MILES 
models (VAZ10) together with those of MAR08 to interpret our 
results. The choice of these models is deliberate. The MILES mod- 
els currently provide one of the best libraries for the study of 
stellar populations at o ptical wav elengths. They are based on the 
MILES stellar library dSanchez-Blazquez et al.l |20 06) which con- 
sists of ~ 1000 stars spanning a large range in atmospheric pa- 
rameters and covering the wavelength range 3525-7500A at 2.3A 
(FWHM) spectral resolution. The spectral resolution, spectral type 
coverage, flux calibration accuracy and number of stars represent a 
substantial improvement over previous libraries used in population 
synthesis models. The MAR08 models on the other hand provide 
colour predictions at infrared wavelengths, and in particular for the 
Spitzer/IRAC 3.6/^m-band data used here. These models differ from 
MILES in that they have incorporated an improved treatment of the 
thermally-pulsing asymptotic giant branch (TP-AGB) phase. The 
main advantage is that both libraries make use of the same set of 
theoretical isochrones JGirardi et alJuOOOf) and therefore allow us 
to generate photometric and spectroscopic predictions on the same 
age and metallicity scales, as well as for the same type of initial 
mass function (IMF). For the purpose of our work we focus on pre- 
dictions for a Kroupa IMF (not corrected for binaries). The spectro- 
scopic measurements have been performed after transforming the 
MILES models to the LIS-14.0A system (see VAZ10). 

Although still somewhat controversial in the details, the im- 
portance o f the TP-AGB ph ase on observed colours has been recog- 
nised (e.g. lMarastonll2005r) . In order to understand the impact the 
TP-AGB phase has on the optical and infrared colours used here, 
we show in Figure lATI the photometric predictions of both sets of 
models up to the reddest band predicted by the MILES models (K- 
band). Given the small difference observed between the K- and 
3.6/^m colours, any conclusion based on A'-band can be extended 
to the Spitzer/IRAC 3.6pm data. The top row shows the V — K 
colour as a function of age for a range of metallicities. The differ- 
ence between the MILES and MAR08 models, at large metallici- 
ties in particular, is significant. While the MILES models predict a 
strong correlation between colour and age, this is essentially non- 
existent for ages above 1 Gyr in the MAR08 models (which show a 
constant behaviour). This result contrasts with the good agreement, 
albeit with a small offset, between the two sets of models for the 
B — V colour (as shown in the middle panel) over the range of 
ages and metallicities displayed by our data. This is confirmed in 
the bottom panel, showing that the predicted l/-band magnitudes in 
both set of models are consistent except for an offset. 

The large similarities between the MILES and MAR08 models 
at optical wavelengths allow us to safely combine the spectroscopic 
prediction from MILES and the photometric ones from MAR08 (as 
shown in Fig. [7}. 

As shown in i |6.1| our sample galaxies display a very tight cor- 
relation between the (V — [3.6]) e colour and a c . Transforming that 
relation into a mass-metallicity relation requires to derive a rela- 
tion between the (V — [3.6]) and metallicity ([M/H]) first. This 
is shown in Fig. lA2l based on the MAR08 models for the range of 
colours and H/3 values exhibited by our sample ((V— [3.6]) e > 2.8 
magandH/3 o <3.0A). 




1.0 

Age (Gyr) 

Figure Al. Compa rison of the c olour and luminosity predictions of the 
iMarigo et al] fe008l> and MILES IVazdekis et ai]|2010l) models as a func- 
tion of age (in Gyr) and metallicity (increasing from —2.32 to +0.22, from 
bottom to top of the figures, with increasing colour or luminosity), for a 
Kroupa initial mass function. 




-0.5 
[M/H] 

Figure A2. (V — [3.6]) vs metallici ty relation from the single stellar popu- 
lation models of lMarigo et al.M 2008). Dashed lines mark the predictions for 
different ages and metallicities. The thick solid line shows the best fit for the 
range of colours and H/3 values defined by our sample (top-right quadrant 
delimited by the horizontal and vertical dashed lines). The fitted relation is 
also indicated. 



APPENDIX B: STELLAR MASS-TO-LIGHT RATIOS 

The effective stellar mass-to-light ratio (7.) is one of the most 
sought after quantities when it comes to determine the stellar mass 
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content or to derive the fraction of dark matter present in a galaxy 
(e.g. after comparison with dynamical esti mates). Pioneering wor k 
relied on optical and infrared colours (e.g. iBell & de Jondl200lh . 
This has been extended over the yea rs to provide new calibra- 
tions based on dif ferent sets of filters JLonghetti & Saraccoll2009l ; 
IZibetti et al J2009Q . The use of full sp ectral fitting methods has also 
been extensive (e.g. lPanter et al.l2007l) , but as a drawback it requires 
some prescription for the amount of dust present in the fitted galax- 
ies. An alternati ve not much explored is the use of line-str ength 
indices (but see lAUanson et alj|2009l ; braves & Faberj|2010l) . The 
advantage is that line-strength s are largely insensitive to the pres- 
ence of dust dMacArthur]l2005l) and onl y a small number of in dices 
is required to get accurate 7, estimates JGallazzi & Belll2009n . 

Ideally one would like to use 7* from models which include a 
proper treatment of the TP-AGB phase. The MAR08 models how- 
ever do not provide such a quantity. It is possible though to estimate 
7» based on the MILES predictions corrected for the colour differ- 
ences of the two model libraries (given their consistency in the opti- 
cal range). For the V-band, for instance, we scale the MILES 7* for 
the small difference between the models' V-band predictions (bot- 
tom panel in Fig. lAll l. Once this is done, it is possible to determine 
7* at 3.6/xm using the MAR08 (V - [3.6]) c colour. 

Detailed studies on the star formation histories (SFHs) of 
early-type galaxies show a wide range of scenarios. It does seem 
however tha t the SFH of a gal axy is largely dependent on its 
stellar mass JThomas et al J 1 20051) . with massive galaxies showing 
single stellar populations while lower mass galaxies display more 
extended star formation. The sample considered here consists 
of early-type galaxies as well as Sa disc galaxies. A priori the 
expected range of SFHs is quite large given the range of single 
stellar population quantities exhibited (see Papers XI and XVII). 
While full, non-parametric spectral fitting is desirable in principle, 
it is not clear whether and how the short SAURON spectral range 
will provide more constraints than those from the H/3, H/3 , 
Fe5015 and Mgfe line-strength indices (i.e. the main features in our 
wavelength range). 

We choose here to use as baseline SFHs different combinations 
of two stellar populations. While a library with different combina- 
tions might not accurately represent the SFH of any given galaxy, it 
will at least show the maximum range of 7, allowed for a given pair 
of H/3 and [MgFe50]' indices. We thus built a library of ~ 5 x 10 6 
models by combining two single stellar populations in the age range 
from 0.1 to 17.28 Gyr (sampled logarithmically) and [M/ H] =[- 
2.32,0.22] (sampled uniformly) for different mass fractions of the 
young component from (i.e. an SSP) to 0.25. We set no restric- 
tion on the metallicities of the two components, and simply impose 
the age of the young one to be lower than the old one. Figure IBTI 
shows the median 7*, for the V- and 3.6/^m-bands, in the H/3 
vs [MgFe50]' space. It also shows the number of models per bin. 
On average our library of models contains at least 1000 predictions 
per bin. Figur JBTI shows that 7, variations depend mostly on H/3 
and that this dependency is stronger in the 3.6^m-band than in the 
V-band. 

We d etermine the best effect ive 7* via a Bayesian approach 
following Kauffmann et al. (2003) (see Appendix A in that paper 
for details). This method results into a probability distribution for a 
given set of plausible models. We take as effective 7* the weighted 
mean value of the distribution and the 68% confidence interval as 
its uncertainty. Measured values are listed in Tables |D3l and ID4I 



£4 
X 



£4 
X 



£4 
X 







12 3 4 

[MgFe50]' (A) 

Figure Bl. Stellar mass-to-light ratios (7*) predictions for the set of mod- 
els produced by the combination of two single stellar populations. Top and 
middle panel show the median M/L ratio in the V— and 3.6/im-band. The 
bottom panel shows a two-dimensional diagram with the number of models 
(in logarithmic scale) falling in a given bin in the H/3 vs [MgFe50]' space. 
For reference, dashed lines mark locations of equal age (along the vertical 
direction) and metallicity (along the abcissae) for a single stellar population. 



APPENDIX C: DEPENDENCE OF KINEMATIC 
SUBSTRUCTURE AND ENVIRONMENT IN SCALING 
RELATIONS 

In this Appendix we show the main scaling relations in the V- and 
3.6/iin-band, colour-coded according to kinematic substructure and 
environment. They are shown here, rather than in the main text, for 
completeness as they do not reveal any particular trend. 



APPENDIX D: PHOTOMETRIC & SAURON QUANTITIES 

Here we present tables with the aperture photometry derived for 
the SAURON sample galaxies in both the V- and 3.6/im-bands, as 
described in the text (see §[4}. We also give tables with the spec- 
troscopic quantities derived from the SAURON data, as described in 
§[5] and Appendix [B] 
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Figure CI. Scaling relations in the V- and 3.6/^m-band for the SAURON sample. Galaxies are colour-coded according to their level of kinematic substructure: 
(CLV) central low velocity, (KDC) kinematically decoupled core, (MC) multiple component, (SC) single component. Elliptical and lenticular galaxies are 
marked with circles, while Sa galaxies are denoted with diamonds. Best-fit relations as shown in J6]and S|2] 
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Figure C2. Scaling relations in the V- and 3.6/^m-band for the SAURON sample. Galaxies are colour-coded according to their environment as indicated. 
Elliptical and lenticular galaxies are marked with circles, while Sa galaxies are denoted with diamonds. Best-fit relations as shown in ij6]and 
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Table Dl. Photometric quantities for the SAURON sample of E/SO galaxies. 



Galaxy 


Env. 


D 


Ref. 


R c ,v 


(Mo,V> 


My 


A V 


R e,[3.6] 


(Me,[3.6]> 


-^"[3.6] 


A[3.6] 


(V - [3.6])e 


ny 


n [3.6] 






(Mpc) 




(arcsec) 


(mag arcsec -2 ) 


(mag) 


(mag) 


(arcsec) 


(mag arcsec -2 ) 


(mag) 


(mag) 


(mag) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


0474 


F 


29.51 


1 


28.08 ±1.65 


20.72 ± 0.09 


-20.87 ±0.12 


0.11 


20.07 ±0.37 


20.72 ±0.11 


-23.94 ±0.11 


0.01 


3.26 ±0.05 


4.0 


4.0 


0524 


F 


23.99 


3 


35.41 ±1.02 


19.90 ±0.07 


-21.75 ±0.20 


0.27 


28.86 ±0.26 


19.90 ±0.08 


-25.04 ±0.20 


0.02 


3.43 ±0.03 


2.8 


2.8 


0821 


F 


24.10 


3 


31.33 ±1.49 


20.09 ±0.08 


-21.30 ±0.17 


0.36 


22.88 ±0.46 


20.09 ±0.10 


-24.31 ±0.17 


0.02 


3.19 ±0.04 


4.0 


1.0 


1023 


F 


11.43 


3 


48.75 ±1.46 


19.65 ±0.05 


-21.08 ±0.16 


0.20 


26.95 ±0.19 


19.65 ±0.08 


-24.08 ±0.16 


0.01 


3.35 ±0.02 


1.0 


3.0 


2549 


F 


12.65 


3 


13.92 ±0.43 


18.84 ±0.16 


-19.38 ±0.28 


0.21 


11.49 ±0.08 


18.84 ±0.19 


-22.56 ±0.28 


0.01 


3.30 ±0.02 


2.9 


2.9 


2685 


F 


16.60 





23.68 ±1.16 


20.02 ±0.10 


-19.95 ±0.49 


0.20 


16.27 ±0.41 


20.02 ±0.14 


-22.94 ±0.49 


0.01 


3.21 ±0.04 


1.0 


3.5 


2695 


F 


32.36 


3 


18.74 ±1.26 


20.05 ±0.12 


-20.86 ±0.40 


0.06 


11.72 ±0.13 


20.05 ±0.19 


-23.96 ± 0.40 


0.00 


3.37 ±0.04 


4.0 


3.1 


2699 


F 


26.92 


3 


12.14 ±0.95 


19.59 ±0.18 


-19.97 ±0.27 


0.07 


7.83 ± 0.08 


19.59 ±0.28 


-22.99 ±0.27 


0.00 


3.29 ±0.05 


1.0 


2.9 


2768 


F 


22.39 


3 


67.97 ±3.11 


20.94 ±0.05 


-21.97 ±0.24 


0.14 


45.49 ± 0.28 


20.94 ± 0.05 


-25.03 ±0.24 


0.01 


3.29 ±0.04 


1.0 


3.5 


2974 


F 


21.48 


3 


28.35 ±1.61 


19.87 ±0.09 


-21.05 ±0.24 


0.18 


20.33 ±0.11 


19.87 ±0.11 


-24.27 ±0.24 


0.01 


3.42 ±0.04 


1.0 


3.4 


3032 


F 


21.98 


3 


19.30 ±1.07 


20.41 ±0.12 


-19.73 ±0.28 


0.06 


10.61 ±0.09 


20.41 ± 0.20 


-22.49 ±0.28 


0.00 


3.08 ±0.04 


3.6 


3.6 


3156 


F 


22.39 


3 


14.78 ±0.44 


20.02 ±0.15 


-19.58 ±0.14 


0.11 


13.40 ±0.14 


20.02 ±0.16 


-22.39 ±0.14 


0.01 


2.89 ±0.03 


2.1 


2.1 


3377 


C 


11.22 


3 


38.28 ±1.99 


20.02 ±0.07 


-20.14 ±0.10 


0.11 


21.80 ±0.22 


20.02 ±0.10 


-22.92 ±0.09 


0.01 


3.08 ±0.04 


1.0 


1.0 


3379 


C 


10.57 


3 


44.89 ± 1.47 


19.47 ±0.05 


-20.90 ±0.11 


0.08 


33.53 ± 0.40 


19.47 ±0.07 


-24.04 ±0.11 


0.00 


3.31 ±0.03 


4.0 


1.0 


3384 


C 


11.59 


3 


28.50 ±0.91 


19.15 ±0.08 


-20.44 ±0.14 


0.09 


20.41 ±0.09 


19.15 ±0.11 


-23.57 ±0.14 


0.00 


3.28 ±0.02 


4.0 


1.0 


3414 


F 


25.23 


3 


32.02 ±1.45 


20.31 ±0.08 


-21.22 ±0.33 


0.08 


20.47 ±0.23 


20.31 ±0.11 


-24.29 ± 0.33 


0.00 


3.31 ±0.03 


1.0 


1.0 


3489 


C 


12.08 


3 


21.51 ±0.81 


18.82 ±0.10 


-20.25 ±0.15 


0.05 


13.38 ±0.14 


18.82 ±0.16 


-23.14 ±0.15 


0.00 


3.15 ±0.02 


1.0 


3.0 


3608 


F 


22.91 


3 


33.64 ±1.52 


20.34 ±0.07 


-21.09 ±0.14 


0.07 


23.98 ±0.15 


20.34 ±0.09 


-24.14 ±0.14 


0.00 


3.24 ±0.04 


4.0 


1.0 


4150 


C 


13.74 


3 


15.88 ±0.52 


19.62 ±0.14 


-19.07 ±0.24 


0.06 


10.10 ±0.09 


19.62 ±0.21 


-22.04 ± 0.24 


0.00 


3.24 ±0.02 


2.6 


2.6 


4262 


C 


15.42 


2 


10.60 ±0.43 


18.53 ±0.21 


-19.54 ±0.06 


0.11 


7.64 ±0.07 


18.53 ±0.28 


-22.71 ±0.06 


0.01 


3.34 ±0.02 


1.0 


1.0 


4270 


C 


16.53 


5 


13.68 ±0.47 


19.74 ±0.16 


-19.02 ±0.09 


0.07 


13.34 ±0.14 


19.74 ±0.16 


-22.21 ±0.08 


0.00 


3.21 ±0.03 


2.2 


2.2 


4278 


C 


16.07 


3 


30.57 ±1.07 


19.48 ±0.07 


-20.97 ±0.20 


0.09 


22.51 ±0.39 


19.48 ±0.10 


-24.14 ±0.20 


0.00 


3.32 ±0.02 


4.0 


1.0 


4374 


C 


18.45 


2 


70.08 ±2.95 


20.17 ±0.04 


-22.39 ±0.06 


0.13 


47.95 ±0.60 


20.17 ±0.05 


-25.46 ±0.05 


0.01 


3.28 ±0.03 


4.0 


1.0 


4382 


C 


17.86 


2 


94.40 ±3.17 


20.58 ±0.04 


-22.55 ±0.06 


0.10 


74.77 ±0.70 


20.58 ± 0.03 


-25.60 ±0.05 


0.01 


3.21 ±0.04 


1.0 


1.0 


4387 


C 


17.95 


2 


11.04 ±0.22 


19.33 ±0.20 


-19.15 ±0.06 


0.11 


10.43 ± 0.08 


19.33 ±0.21 


-22.28 ±0.06 


0.01 


3.18 ±0.02 


2.1 


2.1 


4458 


C 


16.37 


2 


19.94 ±0.96 


20.46 ±0.12 


-19.11 ±0.06 


0.08 


15.21 ±0.14 


20.46 ±0.14 


-21.99 ±0.05 


0.00 


3.07 ±0.04 


3.1 


2.9 


4459 


C 


16.07 


2 


40.97 ±1.48 


20.09 ± 0.06 


-21.00 ±0.07 


0.15 


26.54 ±0.36 


20.09 ± 0.08 


-24.12 ±0.06 


0.01 


3.35 ±0.03 


4.0 


1.0 


4473 


C 


15.28 


2 


26.76 ±0.89 


19.29 ±0.08 


-20.76 ±0.06 


0.09 


22.26 ±0.31 


19.29 ±0.10 


-23.97 ±0.06 


0.00 


3.30 ±0.02 


4.0 


1.1 


4477 


C 


16.53 


5 


46.48 ±2.22 


20.48 ± 0.06 


-20.94 ±0.09 


0.10 


35.87 ±0.69 


20.48 ± 0.06 


-24.10 ±0.08 


0.01 


3.32 ±0.04 


1.0 


1.0 


4486 


C 


17.22 


2 


106.18 ±4.23 


20.60 ±0.05 


-22.70 ±0.08 


0.07 


80.69 ±0.71 


20.60 ± 0.03 


-25.98 ±0.07 


0.00 


3.44 ± 0.04 


1.0 


1.0 


4526 


C 


16.90 


3 


35.70 ±1.29 


19.44 ±0.07 


-21.46 ±0.20 


0.07 


24.03 ±0.22 


19.44 ±0.09 


-24.75 ± 0.20 


0.00 


3.53 ±0.03 


3.3 


3.2 


4546 


C 


14.06 


3 


21.98 ±0.94 


19.03 ±0.10 


-20.41 ±0.20 


0.11 


13.60 ±0.10 


19.03 ±0.16 


-23.48 ±0.20 


0.01 


3.34 ±0.03 


1.0 


3.1 


4550 


C 


15.49 


2 


11.55 ±0.27 


19.00 ±0.19 


-19.26 ±0.05 


0.13 


10.13 ±0.06 


19.00 ±0.21 


-22.30 ±0.05 


0.01 


3.15 ±0.02 


1.8 


1.8 


4552 


C 


15.85 


2 


33.89 ±1.09 


19.43 ±0.07 


-21.21 ±0.06 


0.13 


27.88 ±0.21 


19.43 ± 0.08 


-24.53 ± 0.06 


0.01 


3.43 ±0.02 


4.0 


1.0 


4564 


C 


15.85 


2 


19.35 ±1.11 


19.41 ±0.12 


-20.02 ±0.07 


0.11 


12.97 ±0.21 


19.41 ±0.17 


-23.13 ±0.06 


0.01 


3.34 ±0.04 


4.0 


3.4 


4570 


C 


17.06 


2 


12.75 ±0.56 


18.48 ±0.17 


-20.20 ±0.07 


0.07 


9.17 ±0.08 


18.48 ±0.24 


-23.33 ±0.06 


0.00 


3.34 ±0.03 


1.0 


2.1 


4621 


C 


14.93 


2 


46.06 ±2.99 


19.96 ± 0.06 


-21.22 ±0.07 


0.11 


30.38 ±0.39 


19.96 ±0.07 


-24.32 ± 0.06 


0.01 


3.33 ±0.04 


1.0 


1.0 


4660 


C 


15.00 


2 


11.54 ±0.41 


18.44 ±0.19 


-19.75 ±0.05 


0.11 


8.87 ±0.11 


18.44 ±0.25 


-22.89 ±0.05 


0.01 


3.31 ±0.02 


1.0 


1.0 


5198 


F 


48.78 


8 


18.01 ±0.69 


20.09 ±0.12 


-21.62 ±0.03 


0.08 


15.01 ±0.14 


20.09 ±0.14 


-24.80 ±0.01 


0.00 


3.31 ±0.03 


2.5 


2.5 


5308 


F 


29.43 


1 


9.89 ±0.25 


18.72 ±0.22 


-20.60 ±0.22 


0.06 


7.31 ±0.04 


18.72 ±0.30 


-23.73 ±0.22 


0.00 


3.38 ±0.02 


4.0 


1.7 


5813 


F 


32.21 


3 


55.99 ±2.72 


20.94 ±0.06 


-22.33 ±0.19 


0.19 


39.03 ± 0.62 


20.94 ±0.06 


-25.42 ±0.18 


0.01 


3.33 ±0.05 


1.0 


1.0 


5831 


F 


27.16 


3 


29.19 ±1.70 


20.45 ± 0.08 


-21.04 ±0.17 


0.19 


21.82 ±0.27 


20.45 ±0.10 


-24.11 ±0.17 


0.01 


3.23 ±0.04 


1.0 


1.0 


5838 


F 


23.99 





20.61 ±0.87 


19.29 ±0.11 


-21.18 ±0.49 


0.17 


15.25 ±0.12 


19.29 ±0.14 


-24.43 ± 0.49 


0.01 


3.41 ±0.03 


1.0 


3.8 


5845 


F 


25.94 


3 


4.29 ±0.26 


17.37 ±0.51 


-19.86 ±0.21 


0.17 


3.59 ±0.08 


17.37 ±0.61 


-22.98 ±0.21 


0.01 


3.29 ±0.05 


4.0 


1.0 


5846 


F 


24.89 


3 


76.79 ±3.32 


21.03 ±0.05 


-22.37 ±0.21 


0.18 


64.39 ±0.61 


21.03 ±0.04 


-25.62 ±0.20 


0.01 


3.35 ±0.05 


1.0 


1.0 


5982 


F 


40.36 


1 


24.89 ±1.19 


20.03 ± 0.09 


-21.97 ±0.09 


0.06 


19.65 ±0.19 


20.03 ±0.11 


-25.23 ±0.09 


0.00 


3.38 ±0.03 


1.0 


1.0 


7332 


F 


23.01 


3 


9.21 ±0.15 


18.09 ±0.24 


-20.54 ±0.20 


0.12 


8.14 ±0.04 


18.09 ±0.27 


-23.64 ± 0.20 


0.01 


3.20 ±0.01 


2.0 


2.0 


7457 


F 


13.24 


3 


33.25 ±1.28 


20.59 ±0.08 


-19.63 ±0.21 


0.17 


27.37 ±0.19 


20.59 ±0.08 


-22.55 ±0.21 


0.01 


3.08 ± 0.04 


2.5 


2.5 



NOTES: (1) NGC gala xy number. (2) Envir on ment: Field (F) o r Cluster (C). (3) D istance (Mpc). (4) Original reference for the distance estimate (0: redshift 
distance from NE D . l:ICantiello et alj2007l; 2:lMei et alj2007l ; 3: lTonrv et alj200ll; 4: lReindl et alj2005l ; 5: lMei et alj200l but see (hi) in §33\ 6: 
iTerrv et alj2002t 7: ITu11v et alj2008t 8: lwillick et aljl997l) . See £|3.3| for uncertainties in the measurements, (5-12) Photometric quantities: effective radii 
(Re), mean effective surface brightness ((fi c )), absolute magnitude (M) and Galactic extinction in the V- and 3.6/im bands. (13) Effective (V — [3.6]) 
colour. (14-15) Sersic indices used in our growth curve analysis (see [0). 
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Table D2. Photometric quantities for the SAURON sample of Sa galaxies. 



Galaxy 


Env. 


D 


Ref. 


Ro.v 


<Me,V> 


My 


A V 


R e,[3.6] 


(Me,[3.6]> 


M [3 .6] 


A[3.61 


(V - [3.6])e 


n V 


n [3.6] 






(Mpc) 




(arcsec) 


(mag arcsec -2 ) 


(mag) 


(mag) 


(arcsec) 


(mag arcsec -2 ) 


(mag) 


(mag) 


(mag) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


1056 


F 


32.21 


7 


17.48 ±1.24 


20.17 ±0.13 


-20.57 ±0.40 


0.48 


7.14 ±0.32 


20.17 ±0.31 


-23.84 ±0.40 


0.03 


3.70 ± 0.06 


4.0 


4.0 


2273 


F 


31.62 


6 


27.64 ±1.22 


20.53 ± 0.09 


-21.17 ±0.39 


0.23 


17.78 ±0.65 


20.53 ±0.12 


-24.56 ± 0.39 


0.01 


3.64 ±0.04 


2.9 


6.6 


2844 


F 


25.12 





16.91 ± 1.23 


20.72 ±0.14 


-19.42 ±0.49 


0.06 


7.39 ±0.11 


20.72 ± 0.29 


-22.39 ± 0.49 


0.00 


3.41 ±0.05 


4.0 


2.9 


3623 


C 


11.97 


7 


48.21 ± 0.95 


19.73 ±0.05 


-21.07 ±0.35 


0.08 


31.99 ±0.32 


19.73 ±0.07 


-24.14 ±0.35 


0.00 


3.45 ±0.03 


2.1 


1.5 


4220 


F 


18.79 


7 


20.30 ±0.50 


20.08 ±0.11 


-19.82 ±0.36 


0.06 


13.50 ±0.08 


20.08 ±0.16 


-23.10 ±0.36 


0.00 


3.62 ±0.03 


1.6 


2.1 


4235 


C 


16.53 


5 


16.44 ± 0.49 


19.96 ±0.13 


-19.21 ±0.09 


0.06 


12.30 ±0.24 


19.96 ±0.18 


-22.74 ±0.08 


0.00 


3.71 ±0.03 


1.9 


1.9 


4245 


C 


13.18 





35.95 ± 1.63 


20.77 ±0.07 


-19.60 ±0.49 


0.07 


22.14 ±0.23 


20.77 ±0.10 


-22.50 ±0.49 


0.00 


3.24 ±0.04 


3.5 


3.5 


4274 


C 


19.41 


7 


45.43 ±1.26 


20.56 ±0.06 


-21.16 ±0.36 


0.07 


33.70 ±0.16 


20.56 ±0.06 


-24.45 ± 0.36 


0.00 


3.53 ±0.03 


2.3 


1.6 


4293 


C 


16.53 


5 


52.20 ±1.63 


21.09 ±0.06 


-20.59 ±0.09 


0.13 


31.00 ±0.48 


21.09 ±0.07 


-23.51 ±0.08 


0.01 


3.36 ±0.05 


1.8 


2.1 


4314 


C 


17.38 





46.90 ± 2.24 


20.74 ±0.06 


-20.81 ±0.49 


0.08 


32.80 ±0.59 


20.74 ±0.07 


-23.94 ± 0.49 


0.00 


3.31 ±0.03 


1.0 


3.0 


4369 


F 


20.89 





24.73 ±1.27 


20.38 ±0.09 


-20.18 ±0.49 


0.08 


14.67 ±0.15 


20.38 ±0.15 


-23.01 ± 0.49 


0.00 


3.13 ±0.04 


1.0 


2.7 


4383 


C 


16.53 


5 


14.83 ±0.79 


19.77 ±0.15 


-19.18 ±0.09 


0.08 


9.98 ±0.19 


19.77 ±0.22 


-22.30 ±0.08 


0.00 


3.32 ±0.03 


1.0 


1.0 


4405 


C 


16.53 


5 


18.29 ±0.48 


20.51 ±0.12 


-18.89 ±0.09 


0.08 


12.95 ±0.17 


20.51 ±0.17 


-21.83 ±0.08 


0.00 


3.25 ±0.04 


1.2 


1.2 


4425 


C 


16.53 


5 


15.91 ±0.65 


20.03 ±0.14 


-19.07 ±0.09 


0.09 


12.36 ±0.09 


20.03 ±0.18 


-22.00 ±0.08 


0.00 


3.14 ±0.04 


2.0 


2.0 


4596 


C 


16.53 


5 


59.91 ± 3.05 


20.93 ±0.06 


-21.05 ±0.09 


0.07 


38.24 ± 0.45 


20.93 ± 0.06 


-24.02 ± 0.08 


0.00 


3.25 ±0.05 


1.0 


1.0 


4698 


C 


16.53 


5 


50.64 ± 2.40 


20.76 ±0.06 


-20.85 ±0.09 


0.08 


33.74 ±0.22 


20.76 ± 0.06 


-23.80 ±0.08 


0.00 


3.22 ±0.04 


1.0 


3.0 


4772 


C 


16.53 


5 


42.71 ± 2.34 


21.08 ±0.07 


-20.16 ±0.09 


0.09 


25.26 ±0.31 


21.08 ±0.09 


-23.01 ± 0.08 


0.00 


3.14 ±0.05 


1.0 


1.0 


5448 


F 


34.83 


6 


24.97 ±0.87 


20.88 ±0.10 


-20.81 ±0.18 


0.05 


17.09 ±0.26 


20.88 ±0.13 


-24.00 ±0.18 


0.00 


3.48 ±0.04 


1.8 


3.0 


5475 


F 


30.20 





9.42 ± 0.23 


19.54 ±0.23 


-19.73 ±0.49 


0.04 


7.78 ±0.07 


19.54 ±0.28 


-22.95 ± 0.49 


0.00 


3.36 ±0.02 


1.3 


2.1 


5636 


F 


28.84 





22.89 ±0.89 


21.72 ±0.11 


-19.38 ±0.49 


0.11 


16.65 ±0.24 


21.72 ±0.13 


-21.81 ±0.49 


0.01 


2.85 ±0.07 


1.0 


1.0 


5689 


F 


36.31 





12.63 ±0.35 


19.29 ±0.17 


-21.02 ±0.49 


0.12 


7.91 ±0.04 


19.29 ±0.27 


-24.24 ± 0.49 


0.01 


3.56 ±0.03 


2.1 


2.1 


5953 


F 


33.11 





12.84 ±0.72 


19.53 ±0.17 


-20.61 ±0.49 


0.16 


6.88 ±0.14 


19.53 ±0.32 


-24.03 ± 0.49 


0.01 


3.71 ±0.04 


1.0 


1.0 


6501 


F 


47.86 





20.27 ±1.40 


20.37 ±0.12 


-21.56 ±0.49 


0.29 


14.02 ± 0.33 


20.37 ±0.16 


-24.73 ± 0.49 


0.02 


3.35 ±0.05 


1.0 


5.9 


7742 


F 


22.91 





13.93 ±0.36 


19.17 ±0.16 


-20.35 ±0.49 


0.18 


12.05 ±0.07 


19.17 ±0.18 


-23.41 ± 0.49 


0.01 


3.18 ±0.02 


2.1 


2.1 



NOTES: (1) NGC gala xy number. (2) Envir on ment: Field (F) o r Cluster (C). (3) Dis ta nce (Mpc). (4) Orig ina l reference for th e distance estimate (0: redshift 
distance from NE D : l:ICantiello et al.l2007t 2:lMei et alj2007l : 3: lTonrvet alj200ll : 4: lReindl et all2005t 5: lMei et al.l2007l but see (iii) in JO] 6: 
iTerrv et alj2002l : 7: ITu11v et alj2008l : 8: lwillick et al.ll997l) . See 93.3l for uncertainties in the measurements, (5-12) Photometric quantities: effective radii 
(R e ), mean effective surface brightness ((/i c )), absolute magnitude (M) and Galactic extinction in the V- and 3.6^tm bands. (13) Effective (V — [3.6]) 
colour. (14-15) Sersic indices used in our growth curve analysis (see [|4j. 
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Table D3. Spectroscopic quantities for the SAURON sample of E/SO galaxies. 



Galaxy 


SR/FR 


Kinem. 


U e 


H/3 


Wo 


Fe5015 


Mgfc 


7*,V 


7*, [3. 6] 








(km/s) 


(A) 


(A) 


(A) 


(A) 


(M Q /L V 


(M /L O ) [3 . 6] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


0474 


FR 


MC 


142 


1.85 


2.93 


3.84 


2.96 


o ,-9+0.67 
°° -0.59 


0.78±°- 


0524 


FR 


SC 


225 


1.52 


2.58 


4.26 


3.63 


5 36+ ' 58 
°" :>D -0.60 


i.05±S:ig 


0821 


FR 


SC 


182 


1.65 


2.65 


3.69 


3.22 


4 ^4+0.63 
4 - o4 -0.60 


1 02+ ' 12 
l-V^-O .12 


1023 


FR 


SC 


165 


1.57 


2.55 


3.85 


3.55 


c i 7 +0.56 
°- 1 '-0.57 


1 09+° n 
1 - uy -0.11 


2549 


FR 


MC 


110 


2.01 


3.26 


4.14 


3.10 


9 co+0.90 
^■ b8 -0.65 


53+° ' 17 
u - oc >-0.13 


2685 


FR 


SC 


100 


2.05 


3.09 


3.33 


2.59 


9 SH+0- 74 
2-88_ 069 


72+° 18 
u -' z -0.18 


2695 


FR 


MC 


181 


1.36 


2.19 


3.39 


3.44 


5 01+ 028 


-, 91 +0.03 
I ' zl -0.05 


2699 


FR 


MC 


123 


1.78 


2.83 


3.81 


3.11 


o 91 +0.65 
°- 3 -0.57 


0.87+8:1* 


2768 


FR 


SC 


200 


1.67 


2.67 


3.47 


3.12 


4 9s +0.65 


i.02±8:ii 


2974 


FR 


SC 


227 


1.73 


2.77 


4.01 


3.50 


4 33+O.6O 

^•°°-0.56 


0-89+^2 


3032 


FR 


CLV 


90 


3.85 


5.30 


3.19 


1.74 


-, 01 +0.85 
1 - O1 -0.58 


30+ ' 20 


3156 


FR 


SC 


66 


2.93 


4.16 


2.90 


1.55 


1 45+0.70 


o.4i±8:S 


3377 


FR 


SC 


126 


1.95 


3.10 


3.48 


2.74 


9 co+0.72 


0.69i^ 7 


3379 


FR 


SC 


190 


1.54 


2.60 


3.74 


3.53 


4 qn+0- 58 
^■ au -0.59 


1 06+° n 
1 - UD -0.11 


3384 


FR 


MC 


141 


1.87 


2.98 


4.05 


3.17 


o 4 n +0.66 


071+0.13 

u -' i -0.14 


3414 


SR 


KDC 


191 


1.58 


2.47 


3.49 


3.17 


4 74+0.45 


1 14+O.O8 
l-44_ og 


3489 


FR 


MC 


99 


2.53 


3.65 


3.47 


2.21 


1 Q4+0-78 
1,y -0.54 


46+ ' 20 
u -^°-0.15 


3608 


SR 


KDC 


167 


1.71 


2.78 


3.73 


3.24 


4 lf) +0.64 
4 - iU -0.59 


o.9i± ;il 


4150 


FR 


CLV 


77 


2.69 


3.57 


3.27 


2.00 


2 02+ ' 79 
Z ' Ui! -0.60 


o.52i8:? 7 


4262 


FR 


MC 


164 


1.52 


2.45 


3.45 


3.33 


, 04+0. 42 
4-84_ 043 


-, ic+0.07 
1 - lo -0.09 


4270 


FR 


MC 


125 


1.77 


2.72 


3.64 


2.77 


4 10+ - 68 
^• ±u -0.63 


o.98±8:i^ 


4278 


FR 


MC 


217 


1.65 


2.73 


3.62 


3.63 


4 37+0.63 


095+012 
u - ao -0.12 


4374 


SR 


SC 


261 


1.49 


2.48 


3.61 


3.47 


5 06+ - 48 
a - UD -0.48 


1 14+O.O8 

1 - 1 ^-o.io 


4382 


FR 


CLV 


178 


1.99 


3.09 


3.63 


2.80 


2 93+O 71 


67+ 016 
u -°'_0.16 


4387 


FR 


MC 


98 


1.60 


2.51 


3.58 


3.18 


, 7s +0.49 


1 1 9 +0.09 


4458 


SR 


KDC 


83 


1.64 


2.46 


3.17 


2.88 


4 10+0. 36 
4 - 1S -0.37 


1 ic+0.08 
1 - ±o -0.10 


4459 


FR 


MC 


155 


1.92 


3.07 


3.62 


2.91 


2 go+0-70 
z,aa -0.68 


o.68i8:ii 


4473 


FR 


MC 


186 


1.56 


2.54 


3.82 


3.45 


5 11+0.54 
°- 1± -0.57 


i.ioi8:l? 


4477 


FR 


SC 


147 


1.60 


2.61 


3.70 


3.29 


4 71+0-59 


1 05+° 12 
i-uo_ 12 


4486 


SR 


SC 


268 


1.25 


2.21 


3.52 


3.87 


5 n6+ a3 ° 
OiOD -0.30 


1 91+0.03 

1 - zl -0.05 


4526 


FR 


MC 


214 


1.68 


2.70 


3.89 


3.58 


4 61+ 061 

4 - Di -0.56 


o.96+8:ll 


4546 


FR 


MC 


189 


1.60 


2.49 


3.64 


3.39 


5 02+ - 48 


1 iQ+0. 09 
1 ' lo -0.10 


4550 


SR 


SC 


103 


2.01 


3.06 


3.46 


2.55 


2 07+0.73 
z - 3 '-0.69 


73+ - 17 
u -'-:>_0.i7 


4552 


SR 


KDC 


233 


1.55 


2.64 


4.06 


3.79 


5 03+0.60 
o - uo -0.58 


1 00+° 10 
i.uu_rj.i 


4564 


FR 


SC 


150 


1.62 


2.65 


3.84 


3.47 


4 VK+0.62 
4 -'°-0.59 


101+ 011 
± - u± -0.11 


4570 


FR 


MC 


167 


1.51 


2.45 


3.75 


3.43 


5 94+0-45 
°- /4 -0.47 


1 ic+0.08 
i - ±o -0.09 


4621 


FR 


KDC 


200 


1.50 


2.41 


3.70 


3.56 


5 34+0-41 


1 1 7+O.O6 
J-'- 1 ' -0.09 


4660 


FR 


MC 


181 


1.55 


2.59 


3.72 


3.48 


4.89±°;f 7 


1 07+ ' 11 
^'-O.ll 


5198 


SR 


KDC 


173 


1.60 


2.59 


3.75 


3.43 


4-89+iJ 8 


1 n7+°- n 
^U'-o.ii 


5308 


FR 


MC 


201 


1.52 


2.49 


3.83 


3.60 


5 36+ - 50 


1 iq+0.09 

± - ± °-o.io 


5813 


SR 


KDC 


210 


1.54 


2.61 


3.75 


3.47 


4 8 4+0-59 
4iO4 -0.59 


1 05+° 12 
i ' UO -0.11 


5831 


SR 


KDC 


118 


1.87 


2.96 


3.86 


2.87 


3 4Q+0.68 
°-^ u -0.62 


o.76+8:li 


5838 


FR 


MC 


232 


1.66 


2.67 


3.94 


3.58 


4 76+ - 62 
^•'°-0.57 


0.99+8:11 


5845 


FR 


MC 


237 


1.57 


2.65 


4.19 


3.76 


5 (M+0- 61 
o - uo -0.57 


gg+OiO 
u - aa -o.io 


5846 


SR 


SC 


213 


1.35 


2.36 


3.74 


3.60 


K 4S+0- 37 
5-48_ 042 


1 10+0. 05 
l-18_ oog 


5982 


SR 


KDC 


223 


1.63 


2.66 


4.05 


3.38 


4 79+0.62 
*-' M -0.57 


i.ooi8:ll 


7332 


FR 


KDC 


125 


2.10 


3.27 


3.87 


2.75 


9 „+0.76 


nn + - 17 


7457 


FR 


CLV 


75 


2.26 


3.40 


3.60 


2.40 


9 O7+0.75 
z - z '-0.59 


53+ 019 
U ' M -0.15 



NOTES: (1) NGC galaxy number. (2) Slow/Fast Rotator class. (3) Kinematic substructure as defined in Paper XII. (4) Stellar velocity dispersion within R e ,v- 
(5-8) Line-strength indices within R c y measured in the LIS-14A system. (9-10) Stellar mass-to-light ratios (7*) in the V- and 3.6/im bands. 
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Table D4. Spectroscopic quantities for the SAURON sample of Sa galaxies. 



Galaxy 


SR/FR 


Kinem. 


Ce 


H/3 


n/3 


Fe5015 


Mgb 


l*,v 


7*, [3. 6] 








(km/s) 


(A) 


(A) 


(A) 


(A) 


(M /L )y 


(M O /L ) [3 . 6] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1056 


FR 


SC 


83 


2.85 


3.94 


1.76 


1.24 


1 4V+0.46 


60+ - 18 
U - DU -0.14 


2273 


FR 


MC 


121 


2.72 


3.98 


3.13 


2.27 


1 63+ ' 75 
± - DO -0.50 


40+ a20 
u '^ u -0.14 


2844 


FR 


SC 


93 


2.45 


3.49 


2.56 


1.70 


1 96+ ' 75 


o.64+°; 2 2 ? 


3623 


FR 


MC 


111 


1.69 


2.70 


3.70 


3.27 


4 '3O+0-65 


0.98+°;- 


4220 


FR 


SC 


92 


2.61 


3.68 


3.89 


2.40 


2 f)4+ ' 83 


43+ 017 

u "*°-0.13 


4235 


FR 


MC 


165 


1.84 


2.69 


3.67 


2.97 


4 2Q+ ' 66 


1 00+ 013 
1 - uu -0.13 


4245 


FR 


MC 


97 


1.94 


2.95 


3.43 


2.76 


, q1 4-0.73 


o.8i±g:iS 


4274 


FR 


MC 


153 


1.88 


2.89 


3.51 


2.85 


o cq+0.71 
°- JO -0.62 


0.85±° ;iS 


4293 


FR 


SC 


107 


2.44 


3.53 


3.71 


2.30 


9 i 9 +0.79 
z ' lz -0.56 


0.48+°;?° 


4314 


FR 


SC 


118 


2.13 


3.13 


3.37 


2.60 


2 79+ - 74 


69+ 018 
UiD9 -0.18 


4369 


FR 


SC 


59 


3.39 


4.72 


1.99 


1.30 


i 07 +o.60 

L - yJI -0.41 


36+ 021 


4383 


FR 


SC 


51 


3.05 


4.05 


0.40 


0.89 


i ic+0.40 
1 - 10 -0.25 


RQ+0- 25 

u - DO -o.io 


4405 


FR 


SC 


51 


3.13 


4.36 


2.34 


1.53 


, 91 +0.64 
± - zl -0.37 


0.37±°;S 


4425 


FR 


CLV 


78 


1.89 


2.90 


3.62 


2.78 


o trn+0.70 
00 -0.62 


0.83±g;l| 


4596 


FR 


MC 


156 


1.66 


2.62 


3.51 


3.21 


4 49+0-61 


imt il 


4698 


FR 


KDC 


111 


1.67 


2.59 


3.30 


2.90 


4 20+ ' 53 
4 - ZU -0.55 


l-OS+o'ia 


4772 


FR 


SC 


131 


1.54 


2.50 


3.07 


3.19 


4 27+ ' 42 


-, iq+0.09 
J - la -0.11 


5448 


FR 


MC 


128 


2.19 


3.27 


3.23 


2.65 


2 ^n+ - 79 

Z - OU -0.68 


62+ 020 
U - DZ -0.18 


5475 


FR 


SC 


101 


2.17 


3.35 


4.07 


2.85 


2 4Q+0- 88 
z -^ a -0.62 


51+ 017 
U ' 01 -0.13 


5636 


FR 


CLV 


10 


2.49 


3.64 


2.53 


1.75 


i s1 +0.73 


n cro+0.24 
U - S8 -0.20 


5689 


FR 


MC 


166 


1.98 


3.12 


4.02 


2.93 


2 91+ - 74 
z ' al -0.64 


61+ 015 


5953 


FR 


KDC 


100 


2.86 


3.97 


1.69 


1.29 


1 44+0-46 
^•^-0.39 


60+ 018 
U ' DU -0.14 


6501 


FR 


SC 


190 


1.65 


2.65 


4.06 


3.58 


491 +0.61 


1 00+ 011 
liUU -0.11 


7742 


FR 


KDC 


71 


2.81 


4.07 


2.87 


2.04 


, r-o+0.72 
1 - JO -0.49 


40+ 021 



NOTES: (1) NGC galaxy number. (2) Slow/Fast Rotator class. (3) Kinematic substructure as defined in Paper XII. (4) Stellar velocity dispersion within i? e ,v- 
(5-8) Line-strength indices within i? Ci v measured in the LIS-14A system. (9-10) Stellar mass-to-light ratios (7*) in the V- and 3.6/^m bands. 
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